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ABSTRACT 

We have compiled a significantly updated and comprehensive census of massive 
stars in the nearby Cygnus OB2 association by gathering and homogenising data from 
across the literature. The census contains 169 primary OB stars, including 52 O-type 
stars and 3 Wolf-Rayet stars. Spectral types and photometry are used to place the 
stars in a Hertzprung-Russell diagram, which is compared to both non-rotating and 
rotating stellar evolution models, from which stellar masses and ages are calculated. 
The star formation history and mass function of the association are assessed, and both 
are found to be heavily influenced by the evolution of the most massive stars to their 
end states. We find that the mass function of the most massive stars is consistent 
with a ‘universal’ power-law slope of T = 1.3. The age distribution inferred from 
stellar evolutionary models with rotation and the mass function suggest the majority 
of star formation occurred more or less continuously between 1 and 7 Myr ago, in 
agreement with studies of low- and intermediate mass stars in the association. We 
identify a nearby young pulsar and runaway O-type star that may have originated in 
Cyg OB2 and suggest that the association has already seen its first supernova. Finally 
we use the census and mass function to calculate the total mass of the association of 
16500i!°° Mg, at the low end, but consistent with, previous estimates of the total 
mass of Cyg OB2. Despite this Cyg OB2 is still one of the most massive groups of 
young stars known in our Galaxy making it a prime target for studies of star formation 
on the largest scales. 

Key words: Stars: early-type - Stars: massive - open clusters and associations: 
individual: Cygnus OB2 


1 INTRODUCTION 

Cygnus OB2 is the nearest region to the Sun that 
can lay claim to being amongst the m ost massive clus¬ 
ters or assoc i ations in our galaxy (e.g., IKnodlseder 200(3 : 
lHanson 20031 : IWright &; Drake 2009h . The region is known 
to harbou r many tens of O-type stars and hundreds of O B 
stars (e .g.,|Massev fc T hompson 199 ll ; IComeron et al. 20021 : 

iKiminki et al. 20071 1. though estimates of the total mass 
of the association from both high- and low-mass stars 
range o ver an order of magnitude from (2 — 10) x 
10 4 M.t, (IKnodlseder 2000l : Ih anso n 20031 : IPrew et al. 20081 : 
IWright et al. 2010bl l. The association is itself the most ob¬ 
vious example of recent star formation in the massive 
Cygnus X com pl ex, a regi on of ongoing star formation 
(IReipurth fe Schneider 20081') that span s ~200 pc at a dis¬ 


tance of 1.4±0.08 kpc ( Rvgl et al. 2012 '). The size and prox¬ 
imity of Cyg OB2 therefore make it a template for studies 
of distant, massive star clusters and OB associations in the 
Milky Way and in nearby galaxies. 


The high-mass stellar content of Cyg OB2 has been 
well studied over the past half century since the a ssoci- 
ation was first identified by IMiinch fe Morgan ('19531 1 and 
a number of mid- to late-O stars within it were classified 


by iJohnson fe Morgan f jj)54 ). Further members were iden ¬ 
tifie d and classified by iMorgan et al. f 19541 1 , ISchulte f!956ll , 
and lSchulte ( 19581 ). While many of these classifications w ere 
revised over the following decades (e.g., Walbor n 1973l f. it 
wasn’t until the work of lMassev fc Thompson f 199 il l that a 
more extensive spectral survey of the region was conducted, 
resulting in the identification of 42 O-type stars and 26 B- 
type stars. A decade later the numbe r of kn own O B stars 
incre ased again wi th the works of IComeron et al. 1 20021 1 
and Hanson_[2003j) following the photometric study of 
IKnodlseder f 200ol i . Since then the number of very massive 
O-type stars has not increased significantly, with the more 
recent and de eper surveys uncovering predominantl y B- 
type stars (e.g .. Kiminki ct al. 20071: iKobulnickv et al. 20121 : 

IComeron fe Pasauali 20121 ). This apparent dearth of new O- 
type stars discovered in Cyg OB2, despite deeper observa- 
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Figure 1. Mi d-inf rared 12 //m image of the Cygnus X giant molecular cloud compiled from Wide-Field Infrared Survey 
llWright et al. 201 Oaf) data and spanning 4x3 degrees (100 X 75 pc at the distance of Cyg OB2). Known O, B and Wolf Rayet stars are 
shown as red, green and yellow dots respectively (compiled from the literature as outlined in Section 12.11 The white circle shows the 1 
square degree region studied in this work, which is centred on the “trapezium” of O stars kno wn as Cyg OB2 #8, itself marked with a 
white diamond, and as broadly co-spatial with the Cygnus OB2 Chandra Legacy Survey IWright et al. 2014al) . 


tions, particularly in the near-IR, and through spectroscopic 
observing campaigns, may suggest that our census of mas¬ 
sive O-type stars is nearing completion. 

The age of Cyg OB2 has proved to be a complex sub¬ 
ject, due in no small part to the large spatial size of the 
association a nd the la ck of any clearly recognisable limit 
to its extent. lHanson (2003l l compiled a sample of 85 OB 
stars in Cyg OB2, from which the presence of a number 
of mid O-type dwarf stars and high luminosity blue su¬ 
pergiants led her to conclude that the age of Cyg OB2 
could not be much older than a few million years, and set¬ 
tled on an age of 2 Myrs, with a spread of 1 Myr. The 
evolved giant and supergiant stars present were argued by 
lHanson (20031) to be contaminating non-members, brought 
into the association as the search fo r cluster me m bers wa s 
extended over a larger area (e.g., IComeron et al. 2002h . 
However, studies of A-type stars in Cyg OB2 have re¬ 
vealed a population of 5-7 Myr old stars tow ard and to 
the south of the association rtPrew et al. 200§ 1. while X- 


ray studies of low -mass stars in the a ssociation suggest an 
age of 3-5 Myrs (IWright et al. 2010b! '). If the OB stars in 
Cyg OB2 are younger than the low-mass stars it could 
suggest that the massive stars formed after the low-mass 
stars, a feature observed in some other young clusters (e.g., 
IPovich fe Whitney 2Q10l ~) that has interesting implications 
for theories of massive star formation and feedback. 

The initial mass function (IMF) of high-mass stars pro¬ 
vides a powerful diagnostic of the star formation process 
and is also therefore an area of considerable interest, par¬ 
ticularly as to whether the I MF varies between OB asso¬ 
ciations and star clusters (e.g.. iBastian et al. 20ld ). partic¬ 
ularly given recent suggestions that Cyg OB2 has always 
been a low- density association and did not evolve from a 
star cluster. iMassev fc Thompson (199lh found the high- 
mass IMF in Cyg OB2 to be considera bly fl atter than found 
for other regions in the galaxy, though lKiminki et al. (2007) 
found it to be much steeper than the canonical Salpeter 
IMF llSalpeter 1955l j from an expanded sample, while 
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IWright et al. C2010bl ~l found the lower-mass IMF to be in 
agreement with the Salpeter value. 

Work over the past decade has assembled the data 
needed to form a much fuller view of the OB population 
in Cyg OB2, thanks to a combination of new spectroscopic 
surveys of both p rima ry OB stars and thei r b inary com - 
panions (e.g., iKiminki et al. 20071 : iKobulnickv et al. 2012ll . 
improvements in the stellar effective temperature scale that 
consider line-blanketi ng and no n-LTE (local thermodynamic 
equilibrium) effects (Martins et al. 20051 1 , and the availabil¬ 
ity of revised stellar evolution models for massive stars that 
take into acc ount t he important e f fects o f stella r mass-loss 
and rotation (Ekstrom et al. 20121 ). IWeidner fe Vink f2010l l 
found that masses derived from rotating stellar evolution¬ 
ary models were in much better agreement with dynam¬ 
ical masses for stars in binary systems than were masses 
derived from non-rotating stellar models, which may go 
some way towards addressing previous reports of a non¬ 
standard initial mass function (IMF) in Cyg OB2. Fi¬ 
nally, recent high-spati al resolution optical and near-IR 
Galactic P lane su rveys (Drew et al. 2003 : lLucas et al. 20081 : 
iGroot et al. 20091 ') provide the necessary photometry to re- 
so lve close and ble nded stars in the commonly used catalogs 
of iMassev fe Thompson (T99lll and the 2 Micron All Sky 
Survey 12MASS. ISkrutskie et al. 2006j f - both of which suf¬ 
fer from some degree of blending that is occasionally ignored 
by users of the data. 

Our objective for this work is to bring together all avail¬ 
able information on the massive star population in Cyg OB2, 
gathering spectroscopy from across the literature and se¬ 
lecting the best available photometry for these sources. Our 
hope is that this will allow a quantitative analysis of the star 
formation history, initial mass function, and total mass of 
Cyg OB2. In Section 2 we outline the compilation of spec¬ 
troscopic and photometric data from the literature for our 
census of 169 OB s tars, al most doubling the sample of 85 
OB stars studied bv lHanson~ ('20031 '). In Section 3 we use this 
data to calibrate the extinction towards Cyg OB2 and ver¬ 
ify the extinction law for this sightline that is then used to 
derive individual extinctions towards the massive stars. In 
Section 4 we compile a Hertzsprung-Russell (HR) diagram 
for the population of massive stars, from which individual 
stellar masses and ages are derived. In Section 5 we study 
the age distribution of these stars to infer the star forma¬ 
tion history of the association and in Section 6 we study 
their mass function. In Section 7 we use these results to es¬ 
timate the total stellar mass of the association, and finally 
in Section 8 we summarise our conclusions. 


2 COMPILATION OF THE SAMPLE 

In this section we compile the available photometry and 
spectroscopy from the literature. Figure [T] shows the dis¬ 
tribution of known OB stars in Cyg OB2 and the wider 
Cygnus X region centered on the densest part of the as¬ 
sociation. One of the difficulties with compiling a census 
of Cyg OB2 is that as an association and not a centrally- 
condensed cluster it is not immediately obvious where to 
set the limits of such a study. We focus on a 1 square de¬ 
gree circular area centred on Cyg OB2 #8 - the trapezium 
of O stars at RA 20:33:16, Dec +41:18:45 long regarded 


as being at the h eart of the association ( e.g.. lHanson 20031 : 
IVink et al. 20081) . This is equivalent to a radius of ~14 pc 
at the distance of Cyg OB2. This choice has a number of 
advantages, most notably that it includes the clear over¬ 
density that constitutes Cyg OB2, but does not extend too 
far into the ‘fiel d’ populat ion of OB stars that exists across 
Cygnus X (e.g.. IComeron et al. 20081 ). 


2.1 Spectroscopic data 


The census was compiled by searching the literature for all 
spectroscopically classified stars of spectral type B5 or ear¬ 
lier, a range of spectral types that can provide a solid di¬ 
agnosis of the age, mass function and spatial extent of the 
association. This resulted in 167 stars (22 with known bi¬ 
nary companions) from 18 different publications. To limit 
any potential discrepancies arising from using stars clas¬ 
sified by different authors wherever possible we use clas- 
sifications made by compari son with the stellar atlas of 
IWalborn fe Fitzpatrick (T990h . Fortunately this includes the 
majority of known OB stars in Cyg OB2. When a star is 
classified in this way in multiple papers the practice here 
is to adopt the classification made at the highest spectral 
resolution. 

The majority of star s come from the studies of 
IMassev fe Thompson (T99ll . 50 stars, or 26% of both pri¬ 
mary and secondary stars) and iKiminki et al. 120071 . 78 
stars or 41%), which show excellent agreement with each 
other when their classifications overlap. A further 27 
stars (14%) were taken f rom a series of papers by the 
latter team of authors (IKiminki et al. 2008ll 200911 20121 : 
IKobulnickv et al. 2012) 1. The remaining 37 spectral classifi¬ 
cations (19%) used in this work come from 12 different stud¬ 
ies with between 1 and 7 classifications each, including many 
studies devoted to a single particularly interesting star. Clas¬ 
sifications from these focussed studies were favoured where 
possible as they often uncovered previously unknown binary 
companions that changed the classification of the primary 
star. 


Spectral classifications are presented in Ta¬ 
ble [b 2 with references, known binary companions, 
and source nu m bers u nder the systems initiated by 

I Johnson fe Morg an_J4954 ._th e Sch ulte system, e.g., 

Cyg O B2 # 5), Masse y fe Thompson (T99ll . e.g., MT426), 
and IComeron et al. f2002|. e.g., A26)( We note that the 
star listed as ‘All’_by Comeron_et_aL_f2003) is actually 
MT267 from IMassev fc Thompson (T99ll ). and is listed as 
the latter in this work. By comparing spectral classifications 
of the same star performed by different authors we estimate 
the uncertainty in the spectral types presented here to 
be approximately half a subtype for O and early B-type 
(B0-B1) stars and one subtype for later B-type stars. The 
standard deviation in luminosity class is approximately half 
a luminosity class for O-type stars, one luminosity class for 
early B-type stars (B0-B2), and may be higher for later 
B-type stars, but there is insufficient data for a meaningful 
comparison. These uncertainties are used in all following 
calculations to derive uncertainties on other parameters. 
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2.2 Model atmosphere fits 

In this work physical parameters for stars will mainly 
be based on spectral types from the li terature, using 
the effective tempe rature scales of Martins_et_aL_(200E ), 
iTrundle et al. 120071 '). and iHumphrevs fe McElrov (1984 1. 
However, effective temperatures derived from model atmo¬ 
sphere fits should be more reliable, both because they are 
not quantised in spectral subtypes (or half subtypes), but 
also because they take into account other parameters such 
as stellar mass-loss and chemical abundances. 

We searched the literature and identified 11 stars 
with suitab le model a tmospher e fits from the work s 
of iHerrero et al. 12002 ) and iNegueruela et al. (2008l l. 
We did not inclu de the model atmosphere fits from 
[Herrero et al. fl99fl i since these used mass-loss rates and 
line blanketing analyses that have since been superseded 
(Herrero, private communication). We have also excluded 
stars that have since been found to be binaries with similar 
luminosities, as these will have blended spectra that will 
affect the model atmosphere fits. 

The 11 stars with model atmosphere fits are listed in 
Tablc[T]with the fitted effective temperatures and gravities. 
Compared to the effective temperatures and gravities deter¬ 
mined from the spectral types (also listed in Table [TJ there 
is good agreement within the uncertainties. Given that both 
methods employ the same effective temperature scale, this 
is perhaps not surprising. 

For the remainder of this work for the 11 stars listed in 
Table [T] we used the effective temperatures derived from the 
model atmosphere fits and then calculated intrinsic colours 
and bolometric corrections (BC) for these stars as a func¬ 
tion of effective temperature (instead of spectral type), as 
described below. 


2.3 Photometric data 

Photometric data in the optical and near-IR was gathered 
from the literature to complement the spectroscopy. In the 
optical the primary source of p hotometry was the UBV data 
from lMassev fc Thompson (T99ll i. which is commonly used 
for studies of Cyg OB2. Ho wever, some of these photome tric 
measurements are tagged bv lMassev fe Thompson (199111 as 
‘blended’, including 11 stars in our sample. We therefore 
supplemented the UBV data with Ug'r' ph otometry from 
the Ultra-violet Excess Survey (UVEX, iGroot et al. 2009h 
and converted this onto the UBV system (see Appendix [All. 
For the 11 ‘blended’ sources, 7 had sufficient photometry to 
calculate full UBV photometry while 4 had only U and g' 
UVEX photometry (r 1 was saturated) and so only UB pho¬ 
tometry was calculated. In addition 12 sources in our sam - 
ple were not from the original iMassev fe Thompson (T99ll ) 
samp le (predominantly from the near-IR sample gathered by 
IComeron et al. 20021 1 and we therefore sought UVEX pho¬ 
tometry to homogenise our sample. Of these 12 sources, 9 
had sufficient photometry to derive U and B photometry, 
while the remaining 3 were either saturated in the UVEX 
images or the sources were not observed. For all these sources 
the photometric uncertainties used were propagated from 
the UVEX photometric uncertainties and the dispersions on 
the transformations used. 

Near-IR photometry was taken from 2MASS for the ma¬ 


jority of stars. In view of the low spatial resolution of the 
2MASS observations that could lead to blending in some 
cases we sought replacement photometry from the UKIDSS- 
GPS (United Kingdo m Infrared Deep Sky Survey Galactic 
Plane Survey. ILucas et al. 20081 ). but were only able to ob¬ 
tain unsaturated photometry for 3 sources. For these sources 
we replaced their 2MASS photo metr y by UKIDSS photom- 
etry using the transformations in ILucas et al. (20Ci5h 

The final photometric data used in this study is pre¬ 
sented in Tabic [Bl] on the UBV and 2MASS JHK S photo¬ 
metric systems. The vast majority (93% or 155 out of 167) 
of the primary stars have optical photometry in at least two 
of the UBV bands, and all of the targets have near-IR pho¬ 
tometry in J, H, and K s . 


2.4 Binary stars and their influence on stellar 
parameters 


Unresolved binary systems can lead to the miscalculation 
of the properties of the primary star. The contribution of 
light from a binary companion, which in involved systems is 
cooler than the primary, will make the primary appear cooler 
and more luminous, which could potentially lead to the star 
being diagnosed as more evolved and higher-mass. The bi¬ 
nary population of massiv e stars in Cyg OB2 has recently 


nary population ot massive stars m Uyg Uhz nas recently 
been jwell stud ied (e.g., Kiminki et al. 2007ll 200911 20121 ; 


ICaballero-Nieves et al. 2014l f~ and 22 of our 167 primary 

stars have secondaries with known spectral types, although 
it is likely that there still remain a number of undiscovered 
binary companions in this sample. 

For objects with known companions it is possible to 
correct for the light of the secondary when using the ob¬ 
served photometry as long as the spectral type of both 
stars is known. This is relatively simple when the luminos¬ 
ity ratio of the binary system is known in all the relevant 
bands, or if the bolometric luminosity ratio is known and 
the stars are of similar spectral type (e.g., Cyg OB2 #5 
is an 071+061 binary with a bolometric luminosity ratio 
of 3.1, iLinder et al. 20091 1. For other binary systems this is 
trickier and we must make an assumption about the rela¬ 
tive in-band fluxes of the two stars, for which the simplest 
approach and that adopted here is to assume that the ratio 
of fluxes is equal to the ratio o f i ntrinsic fluxes according 
to s pectral type from [Martins et al. 1 2005h for O-type stars 
and IConti et al. ( 2008h for B-type stars. 

Doing this for our sample we find that the typical (rms) 
adjustment to the colour of the primary star is very small, 
0.006 mag, having a negligible effect on the extinctions. 
This is because the strongest effect occurs when two stars 
have significantly different intrinsic colours, but in such sys¬ 
tems the stars usually have very different luminosities and 
therefore the influence of the secondary is small. The rms 
increase in the observed magnitude is more significant at 
0.5 mag, equivalent to a 37% reduction in the bolometric 
luminosity. This luminosity reduction (magnitude increase) 
varies from 8% (0.09) for Cyg OB2 #11 (O5If+B0V) to 50% 
(0.75) for very similar pairs of stars such as Cyg OB2 #73 
(08111+08111) and B17, an 071+091 binary. 

Cyg OB2 does conta in one known triple system, MT429, 
which IKiminki et al. (20~12h find to be a B0V+B3V+B6V 
system. For simplicity we consider only the two brightest 
components of the system because the B6V star is too faint 
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Table 1. Model atmosphere fits used in this work 


Name Spectral types Model atmosphere fit 



Type 

Ref 

Teff(K) 

log g 

Pff (K) 

VTeff ( K ) 

log g 

°logg 

Ref 

#2 

B1I 

MT91 

24,300 

- 

28,000 

1000 

3.21 

0.1 

H02 

#4 

07111 

K07 

36,077 

3.61 

35,500 

1000 

3.52 

0.1 

H02 

#8C 

05111 

K07 

40,307 

3.69 

41,000 

1500 

3.81 

0.1 

H02 

#10 

091 

K07 

31,368 

3.23 

29,000 

1000 

3.11 

0.1 

H02 

#11 

O5If±B0V 

K12a 

38,612 

3.57 

37,000 

1500 

3.61 

+0.15 

-0.10 

H02 

MT267 (All) 

07.5111 

K12a 

35,019 

3.59 

36,000 

1500 

3.6 

0.2 

N08 

A15 

071 

N08 

34,990 

3.40 

35,000 

1000 

3.2 

0.2 

N08 

A24 

06.5111 

N08 

37,134 

3.63 

37,500 

1500 

3.6 

0.2 

N08 

A26 

09.5V 

N08 

31,884 

3.92 

35,000 

1000 

3.9 

0.2 

N08 

A33 

B0.2V 

N08 

29,390 

- 

31,000 

1500 

4.0 

0.2 

N08 

A38 

08V 

N08 

34,877 

3.92 

36,000 

1500 

4.0 

0.2 

N08 


The 11 stars with model atmosphere fits used in this work, with effective temperatures and gravities taken from the fits and calculated 
from the spectral typ es lis ted. Effective tempe ratures and gravities calculated as a function of spectral type using the data from 
[Martins et al. ( 20051) and [Trundle et al. (20071 1 where possible. Uncertainties on the spectral type parameters are estimated from the 
spectral type and luminosity class uncertainties estimated in Section |2. Hand are approximat ely ±1000 K (approximatel y ±3000 K for 
the early B-type stars) and ±0.1 dex. Model atmosp here fits taken fr omlHerrero et al. ( 20021 H02) a ndlNegueruela et al. l200Sl . N08), 
and spectral types from lMassev fe Thompson (1991] . MT911. iKiminki et al. (20071 . K07). and iKobulnickv et £±*720121 . K12a.i. 


to affect the colour or luminosity of the primary star. We 
also note that while Cyg OB2 #22, an 03If±06V binary 
llWalborn 20021 ). has been resolved in high angular resolu¬ 
tion ima ges wit h a separati on of ~ 1.5" between t he two com¬ 
ponents (lMason et al. 2009l : lMafz Apellaniz 20101 ). these ob¬ 
servations do not provide sufficient photometry to place both 
components on the HR diagram separately, therefore we only 
consider the primary star in this work (after applying the 
corrections for the secondary noted above). 


3 THE EXTINCTION TOWARDS 
CYGNUS OB2 

In this section the spectral types and photometry are used to 
re-derive the form of the extinction law towards Cyg OB2 
and then calculate individual extinctions for all the stars 
in our samp le. To achieve th is we use the unreddened 
colours from [Martins fe Plez 120061) for O -type stars (the 
“observational” T e ff scale) and from lFitzgerald (T970I) and 
iKoornneef f 19831 ) for B-type stars, adjusted for B0-B1 stars 
to provide a smooth transition with the O star colours. 


2.5 Completeness 

It is difficult to estimate the completeness of a sample that 
is compiled from so many different sources of spectroscopy, 
each of which is based on a different method of photomet¬ 
ric selection, area of study, and spectroscopic completeness. 
However, one possible indicator of the level of completeness 
of our sample is the fact that the most recent spectroscopic 
studies (e.g., iKiminki et al. 20071 : iKobulnickv et al. 2012i ) 
have discovered only B-type stars and not new O-type stars, 
possibly suggesting that the majority of O-type stars in 
Cyg OB2 have already been detected. The fact that the 
most recent studies have foun d stars wit h similar extinc¬ 
tions to previous studies (IComeron fc Pasauali 20121 ) also 
suggests that a previously undiscovered population of highly 
obscured (Ay > 10 mag) O-type stars does not exist. Based 
upon this we will make the assumption that our census of 
massive stars is complete for O-type stars in Cyg OB2, which 
is broadly equivalent to being complete down to 15 Mq for 
zero-age main-sequence stars, or 20 Mq if one takes into 
account some further evolution. In a future paper we will 
attempt to test this assumption about our completeness us¬ 
ing new spectroscopic observations of Cyg OB2. 


3.1 The extinction law towards Cygnus OB2 


To determine the form of the extinction law towards 
Cyg OB2 we selected a high-confidence subset of 18 
stars from within our sample to which we perform 
spectral energy distribution (SED) fitting. We selected 
only single O- type main -sequence stars classified by 
iMassev fe Thompson 1199ll) and with no evidence of hav¬ 
ing a similar-luminosity binary companion. By using only 
O-type stars we reduce the uncertainty due to the effective 
temperature and the binarity of the stars, both of which 
increase towards later spectral types. In addition we re¬ 
quired all stars to have been detected in, and have low errors 
(< 0.05 in UBV and < 0.02 in JHK„) in, all six photomet¬ 
ric bands. These criteria resulted in the sample of 18 stars 
in Table [2] 

Reddened stellar spectra w ere c a lculated usi ng CM- 
FGEN non-LTE model spec tra dHillier fe Miller 19981) . as 
used bv lMartins et al . (2005!), redde ned using the extinction 
curves presented by Fitzpatrick fe Massa 120071 ) as a func¬ 
tion of Rv and Ay. Reddening coefficients were determined 
using Kitt P eak UBV filt er prof il es (app ropriate for the ob¬ 
servations oflMassey fe Thompson 19911 ) and 2MASS JHK S 
filter profiles llSkrutskie et al. 20061) . Combining these with 
th e intrinsic c ol ours for O-type dwarf stars presented by 
iMartins fe Plez (20061 ) results in reddened photometry as a 
function of spectral type, Rv, and Ay. 
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Table 2. Stellar sample used to determine the extinction law 
towards Cyg OB2. 


Spectral type Number Ry Ay 


06V 

005 

2.95 

± 

0.04 

6.00 

± 

0.04 

09V 

#14 

2.88 

± 

0.04 

4.49 

± 

0.03 

08V 

#6 

2.89 

± 

0.04 

4.48 

± 

0.03 

08.5V 

#17 

2.93 

± 

0.04 

4.88 

± 

0.03 

08V 

376 

2.89 

± 

0.04 

4.85 

± 

0.03 

08V 

390 

2.83 

± 

0.04 

6.49 

± 

0.04 

08V 

455 

2.85 

± 

0.04 

6.01 

± 

0.04 

09.5V 

#23 

2.90 

± 

0.04 

5.06 

± 

0.04 

08.5V 

#8D 

2.86 

± 

0.04 

5.07 

± 

0.04 

07.5V 

#24 

2.94 

± 

0.04 

5.59 

± 

0.03 

08V 

485 

2.87 

± 

0.04 

5.25 

± 

0.03 

08.5V 

507 

2.89 

± 

0.04 

5.34 

± 

0.04 

07.5V 

534 

2.88 

± 

0.04 

6.26 

± 

0.04 

08V 

#74 

2.89 

± 

0.03 

6.36 

± 

0.04 

07V 

611 

2.97 

± 

0.04 

5.51 

± 

0.04 

09V 

#41 

2.84 

± 

0.04 

6.04 

± 

0.04 

09V 

#75 

3.08 

± 

0.04 

5.48 

± 

0.04 

07V 

#29 

3.02 

± 

0.04 

5.35 

± 

0.03 


All stars classified bv lMassev &; Thompson f!99ll ) with source 
numbers from that work, except those noted with #, which are 
labelled on the Schulte system. Note that while MT005 and 
MT611 have been observ ed to have binary companions 
(ICaba ll ero-Nieves et al. 2014f) , the companions are significantly 
fainter than the primary stars (<51/ = 2.8 and SK = 4.9, 
respectively) and are therefore unlikely to affect the SEDs of the 
primary stars. 


X 2 fits were performed between the observed photom¬ 
etry and the reddened model photometry allowing Ry to 
vary in steps of 0.01 in the range 2.5-3.5 and with Ay vary¬ 
ing between 4.0-7.0 mag in steps of 0.01 mag. The distance 
modulus was left as a free parameter to account for varia¬ 
tions in the intrinsic distance and luminosity of stars. The 
best fits are shown in Figure [2] and listed in Table [2] 

The weighted mean extinction law from our \ 2 fits is 
Ry = 2.91 d= 0.06 (taken as such for the remainder of this 
work), in good agree ment with pr evious studies based on 
optical colours (e.g., lHanson 2003i . measured Ry = 3.0). 
Ho wever t his va lue is in disagreement with the findings of 
IComeron fc Pasquali (201211 who use a combination of op¬ 
tical _and near -IR photometry to derive Rv = 4.0 using 
a ICardelli et al. f!989ll extinction law for OB stars over a 
wider area across the Cygnus X molecular cloud. This dis¬ 
agreement is likely explained by there being two different 
extinction laws at work along different sight lines: a ‘stan¬ 
dard’ Ry ~ 3 extinction law towards Cyg OB2 (due to ob¬ 
scuration from the Cygnus R ift, which lies predominantly in 
the foreground, iGuarcello et al. 20121 ). and an ‘anomalous’ 
Rv ~ 4 extinction law away from Cyg OB2, possible due 
to a change in the grain size distribution favouring larger 
dust grains within the Cygnus X molecular cloud (which 
has already been dispersed from the vicinity of Cyg OB2, 
see Figure [l]). A simi lar picture of two extinction laws has 
been observed by IPovich et al. ('201111 towards the Carina 
molecular cloud. 

For reference this new extinction law is presented in 
Table [3] for extinctions between Av = 4.0-7.0 (the typi¬ 
cal range observed in Cyg OB2, see Section 13.21) and com- 



Figure 2. Ry versus Ay for our high-confidence sample of 
18 single O-type dwarfs derived using y 2 fits between the 
observed p hotometry and reddened model spectra using the 
iFitzpatrick fc Massa ('20071 1 reddening curves, lcr error bars are 
shown derived from the y 2 fits. The weighted mean Ry = 
2.91 ± 0.06 is shown as a dashed line. 


Table 3. Extinction law towards Cyg OB2 


Band (A) 


A x /Ay 



This work 


Hanson 

o 

II 

Ay = 5.5 

o 

[ - 

II 

> 

U (3372 A) 

1.593 

1.599 

1.605 

1.600 

B (4404 A) 

1.321 

1.318 

1.316 

1.333 

V (5428 A) 

1.000 

1.000 

1.000 

1.000 

J (1.27 pm) 

0.233 

0.234 

0.235 

0.282 

H (1.67 fi m) 

0.138 

0.139 

0.140 

0.175 

K s (2.16 /rm) 

0.0840 

0.0845 

0.0850 

0.125 


New extinction law for Cyg OB2 calculated using a 
IFitzpatrick fc Massa (’20071 1 Ry = 2.91 extinction curve and 
various extin ctions. The previous Ry = 3.0 extinction law 
presented by lHanson ('20031 ) is given for reference. 

pared to that found previously bv lHanson" 12003 ) from com- 
parison o f UBVJ H K„ c olour excesses compared to the 
iLeieune fc Schaerer (20011) intrinsic colours. The extinction 
law is relatively unchanged in the optical, but has changed 
significantly in the near-IR, due mostly to shifts in the near- 
IR i ntrinsic colours of O-type stars between th e wo rks of 
ILeieune fe Schaerer f200H ) and I Martins fe Plez (20061 ). 

3.2 The extinction of Cygnus OB2 members 

Using the IFitzpatrick fe Massa (2007ll Ry = 2.91 extinc¬ 
tion law derived above we calculated individual extinctions 
for the 164 primary O and B-type stars in our sample (not 
including the WR stars) from the comparison between ob¬ 
served and intrinsic colours. We favoured the B — K 3 colour 
for this due to its long baseline and the availability of the 
relevant photometry for the vast majority (160 or 98%) of 
our sources. For 1 star B-band photometry was unavailable 
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A v 

Figure 3. Histogram of the extinctions calculated for all 164 
massive OB stars in this work with la error bars derived from 
MC simulations. 



Figure 4. Voronoi tessellated extinction map (one star per 
Voronoi cell) for for all 164 massive stars in Cyg OB2 with the 
extinction indicated by the greyscale colour. The circular border 
of our study is approximately at the border of this image. 


so we used the V — K s colour and for 3 stars only near- 
IR photometry was available so we used the J — K s colour. 
Uncertainties in the derived extinction were calculated us¬ 
ing a Monte Carlo (MC) simulation taking into account the 
uncertainties in the photometry, spectral classification, the 
intrinsic colours (assumed to be half the difference between 
the colours of neighbouring subtypes), and the reddening 
law of Rv = 2.91 ± 0.06. The typical uncertainty in the 
extinction is ~0.3 mag in the U-band. Extinctions and un¬ 
certainties are presented in Table iBlI 

Figure [3] shows the extinction distribution, which varies 
from Ay = 2.2 mag for MT140 (an 09.51 star on the western 
side of the association) to Ay = 10.2 mag for the central blue 
hypergiant Cyg OB2 #12 (a B3.5Ie star, an d one^of the most 
luminous known s tars in the Milky Wav. lMorgan et al. 19541 ; 
IClark et al. 20l3) . though the majority lie in the range 
Ay = 4-7 mags with a median of Av = 5.4 mag and 
25% and 75% quartiles of 4.7 and 6.5 mags. Error bars on 
the extinction distribution histogram were derived from the 
individual uncertainties using a MC simulation. The large 
spread in extinction compared to the typical uncertainty in 
Av suggests that the spread is real and is either due to fore¬ 
ground extinction from the Cygnus Rift or intra-association 
extinction within Cyg OB2. 

Figure [|] shows a spatial map of the extinction across 
Cyg OB2 in the form of a Voronoi tessellation. The ex¬ 
tinction distribution varies smoothly across the association, 
from Ay ~ 4-5 mag in the n orth-west ( where a known ‘red¬ 
dening hole’ exists. IReddish et al. 19671 ) to Av ~ 6-7 to the 
south of the association and on the outskirts. 

These results are in good agr eeme nt with previous stud¬ 
ies of both the high- ( Hanson 2003) and low-mass popu¬ 
lations llWright et al. 2010bl h albeit with a smaller spread 
than the latter, likely due to the higher precision afforded 
by having known spectral types and a calibrated extinc¬ 
tion l aw. The extinc tion is larger than that derived by 
[Guarcello et al. (2012i) from their riz study by ~1 mag 
(they derive a median extinction of Ay = 4.3 mag), 
though we observe the same spatial variation in extinc¬ 


tion that they d o. The an omalously low extinction derived 
by iGuarcello et al. (20121 ), particularly compared to many 
other studies, may be due to the difficulties deriving extinc¬ 
tion from isochrone fits in the riz filter system that requires 
a complex transformation from BVI colours. 

The high extinction for Cyg OB2 #12, > 2 mags larger 
than any of the stars surrounding it, is particularly appar¬ 
ent. Some of this extinction may be due to circumstellar ma¬ 
terial given the advanced evolutionary stage of this object 
and that material is known to surround other very massive 
stars (e.g.. IWright et al. 2014c| ), however ldark et al. (2012i ) 
could find no evidence for a near- to mid-IR excess in the 
spectral energy distribution of this object, arguing against 
a significant amount of warm circumstellar dust. It is also 
possible that the SED we have assumed is wrong, given the 
unusual character of the star. 


4 HERTZPRUNG-RUSSELL DIAGRAM FOR 
CYG OB2 

In this section we place the stars on the HR diagram and 
compare their positions with both rotating and non-rotating 
stellar evolutionary tracks to derive stellar masses and ages. 
Cyg OB2 also includes 4 evolved massive stars (such as Wolf 
Rayet stars) whose positions on the HR diagram and result¬ 
ing stellar properties cannot be derived in the same way as 
the less evolved OB stars. We have gathered the properties 
of these stars from the literature, adjusting for a distance of 
1.45 kpc if necessary, and derive additional physical proper¬ 
ties (such as upper or lower limits on stellar ages of masses) 
by comparison with the evolutionary models used in this 
work or from our current, but far from complete, under¬ 
standing of the evolution of massive stars. These stars and 
their properties are included in Table [Bll and included in all 
later analysis where possible. 
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Figure 5. HR diagram for Cyg OB2 showing all the stars compiled in this work colour-coded by their luminosity class (dwarfs are black, 
subgiants are blue, giants are green, bright giants are orange and supergiants are red, while WR stars are shown with a star symbol) and 
with source num bers n oted according to the scheme adopted in this paper. The figure shows isochrones and evolutionary tracks from 
iLeieune K: Schaerer 120011 ) with ages and masses labelled. The effective temperature scales used for main-sequence O- and B-type stars 
are shown along the top of the diagram and mean la errors bars for each subtype are shown to the bottom left. 


4.1 Calculation of Lbol and T e // 


Bolometric luminosities were calculated using the observed 
K s band magnitudes because of the high extinction to¬ 
wards Cyg OB2 and the typical uncertainty of 0.3 mag 
in Ay resulting from this. The considerably lower un¬ 
certainty in Ak compared to Ay is preferable despite 
the higher uncertainty in BCk over BCy (due to the 
strong sensitivity of the spectral ener gy dist rib ution to win d 
parameters at longer waveleng t hs, Martins fe Plez 20061) . 
BC were tak en from iMartins fe Plez (20061 ) for O-type 


stars, 


Crowther = et i j ; L_(200^) for B-type supergiants, and 


[Humphreys fe McElrov (19841) for less luminous B-type 
stars, although there is very little variation with lu¬ 
minosity class. Effective temperatures were determined 
as a function of spectral type using the tabulations of 


Maxtins_e£_aL_[2005, “observed” scale) for O-type stars, 
Tnmdle_el L jr]^(^007 ) for ear ly (BO-3) B-type stars and 
Humphreys fe McElrov (19841) for later B-type stars, all of 
which show a good agreement between each other. 

The accepted distance to Cyg OB2 has been re¬ 
duced over the last few deca des due to improved obser¬ 
vations and calibration data. [Reddish et al. (19661 ') origi- 
nally placed the association at a distance of 2.1 kpc, while 
iMassev fe Thompson f!99ll ) estimated a nearer distance of 
~1.75 kpc (DM=11.20) following their spectroscopic study 
of massive stars in the association. The r ecent recalibra - 
tion of OB star absolute magnitudes led [Hanson (20031 ) 
to determine an average distance of 1.2 kpc (DM=10.44) 
but concluded based on inconsistencies with current stellar 
wind theory that this was too low and suggested a value 
of 1.45 kpc (DM=10.80). This is the value we use in this 
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work, which is in agreement with the p arallax distance o f 
1.40±0.08 kpc for the Cygnus X complex (iRygl et al. 201 jf . 

The uncertainties in our values of T e // and Lbol were 
derived from MC simulations. For t he former we use th e 
quoted u ncertai nties in T e f f from MarMns_etjiLJ^005), 
iTrundle et al. 120071 1 , and iHumphrevs fc McElrov 11984 ), 
while for the latter we take into account the uncer¬ 
tainties in the observed photometry, derived extinctions, 
BCs (a = 0 .1 mag for O-type stars in the A'-band, 
iMartins fc Plez 20061 . and assumed to be a = 0.2 mag for 
B-type stars), and spectral classification. For the distance to 
Cyg OB2 we assume that the distance of 1.45 kpc is correct 
(small variations in this value up to a few hundred parsecs 
do not significantly affect the results of this paper, though 
larger values do). 


4.2 The Hertzsprung-Russell diagram 

The new HR diagram for Cyg OB2 showing 167 stars 
is shown in Figure [5] compared to the isochrones and 
evolution a ry tracks from the non-rotating models of 
iLeieune fc Schae rer 120 0111 and in Fig ure 151 compared to the 
rotating models of lEkstrom et al. 120121 '). The non-rotating 
models at a metallicity of Z = 0.02 are shown partly 
for comparison with prev ious studies (e.g., IHanson 2003 : 
IComeron fc Pasauali 2012 1. and also to contrast with the 
more recent, rotating models at the revised solar metallicity 
of Z = 0.014. An indication of the uncertainties on Lboi and 
T e ff are shown in Figures [5] and [6] for each spectral sub¬ 
class (where possible) and in Figure [6] for individual stars to 
illustrate the spread they induce in the HR diagram. For O- 
type stars the uncertainties are typically small, ~0.02 dex in 
log Teff and ~0.1 dex in log L/Lq, though for the B-type 
stars the errors rise significantly to ^0.07 dex in log T e // 
and ~0.2 dex in log L/Lq due to classification uncertain¬ 
ties and the larger difference in log T e ff between spectral 
subtypes. 

As noted by previous a uthors (e.g., 
[Massey fe Thompson 199H ; IHanson 2003lf there is a 
well-populated main sequence extending up to 05. 5\E for 
which there is now the data to extend it down to mid-B 
spectral types. 


4.3 A large luminosity spread amongst the B stars 

Many evolved massive stars are evident, both at masses M > 
25 Mq where such evolution can occur within 2-10 Myr, 
and at lower masses where evolved B-type stars occupy 
positions in the HR diagram suggesting ages > 10 Myrs. 
IHanson 12003 1 noted in her sample a large luminosity spread 
in the Bl-2 dwarfs and this is also apparent here, although 
it is now more concentrated towards the main-sequence 
(most likely because the massive stars discovered recently 
have been faint er as t he depth of s pectroscopic observations 
has increased). IHanson (2003 1 and INegueruela et al. (20081 1 
have suggested that their position in the HR diagram may 


1 The 05V star A37 falls outside of our survey area and therefore 
is not included in this study, and may not a ctually be a member 
of Cyg OB2, as suggested by IHanson (200j) . 


be due to them being foreground contaminants from a sep¬ 
arate population. However we note that for their luminosi¬ 
ties to be consistent with foreground, main-sequence B-type 
stars would require some of the most over-luminous stars to 
lie at distances as low as ~500 pc, which would be inconsis¬ 
tent with their high extinctions (the five most over-luminous 
B-type stars all have Ay > 5.5 mag) and the estimated 
distance to t h e extin guishing Cygnus rift of 800-1400 pc 
dGuarcello et al. 20121 : Sale et al. 20141 ). 

Alternatively the apparently large luminosity spread 
amongst the B-type stars could be due to the large uncer¬ 
tainties in Lboi and T e ff at these spectral types that are 
smearing out the main sequence in this range. Observational 
uncertainties in spectral type combined with the steep re¬ 
lationship between T e // and spectral subtype in this range 
creates large uncertainties in T e ff amongst the early B-type 
stars. Comparing the observed spread in the HR diagram 
with our estimated uncertainties suggests there is no signif¬ 
icant evidence favouring a foreground population of B-type 
stars or a co-spatial population of mature (>10 Myr) stars 
to explain the observed spread. This argument is supported 
by the fact that the seemingly high-luminosity B-type stars 
have a similar spatial distribution to the low-luminosity B- 
type stars, suggesting they are unlikely to be a separate 
population. 


4.4 Possible foreground and background 
contaminants 

A small number of stars fall significantly below the main se¬ 
quence in the HR diagram such as MT170 (B5V, log L/Lq = 
1.74), MT426 (BOV, log L/Lq = 2.51), and MT427 (B4III, 
log L/Lq = 2.19). MT170 and MT427 appear to be back¬ 
ground sources, most likely at distances of ~ 2.19-3.5 kpc 
if they are on the main-sequence (a number of Hu regions 
and massive young stell ar objects have been identified at 
similar distances bv lXu et al. 20130 . Their moderate extinc¬ 
tions {Ay = 3.6 and 3.8 mag, respectively) are acceptable 
given that they are projected against the ‘ reddening hole’ in 
the north-west of the association (see also iDrew et al. 20081 . 
who find evidence of a population for A-type stars at such 
distances in this direction). 

The low luminosity of MT426 is harder to explain 
as it would imply that the star is ~10 kpc distant yet 
has an extinction of only Ay ~ 4 whilst being projected 
against the centre of the association where the extinc¬ 
tion is typicall y Ay = 6-7 mag ( Figure [4]). MT426 is 
listed by iMassev fe Thompson (T99ll i as having ‘blended’ 
photometry (see Appendix 0 with V = 14.05 mag and 
B — V = 1.95 mag. These values have led previous authors 
to derive reasonable extinctions and luminosities for this 
star, yet such values would imply a AVband magnitude of 
~8.4, yet the source is not detected by 2MASS. Inspection 
of the UVEX and UKIDSS images reveals a second point 
sour ce that is highly blended in the lower spatial resolu¬ 
tion IMassev fe Thompson (T99ll i and 2MASS images with 
its neighbour MT425. The UVEX and UKIDSS photome¬ 
try, which we have used, does not appear discrepant and 
we believe it to be sound. It is possible that the spectral 
type of B O V is erroneous , given the blending reported by 
IMassev fe Thompson 119911 1 and the fact that those authors 
report an identical spectral type for its neighbour (though 
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Figure 6. HR diagram for Cyg OB2 showing all the stars compiled in this work as per Fig ure [5l with the difference that the main panel 
shows isochrones and evolutionary tracks for rotating {vi n i/v cr it = 0.4) models from lEkstrom et al. (20121 , smoothed for clarity) with 
the ages and masses labelled. The bottom left figure shows the error bars for each source and the isochrones and evolutionary tracks. 
The bottom right figure shows the distribution of all stars without labels and evolutionary tracks. 


© 0000 RAS, MNRAS 000, 000-000 















































































The Massive Star Population of Cygnus OB2 11 


iKiminki et al. 20071 . report the same spectral type, although 
it may also have suffered blending). 

Due to the unresolved nature of MT426 and the fact 
that MT170 and MT427 appear to be background sources 
we argue that all three are non-members of Cyg OB2 and ex¬ 
clude them from further analysis of the association. They are 
noted in Table iBlI as ‘Unresolved’ or ‘Background’ sources, 
and physical properties are not derived for them. 

Finally we note that the membership status of 
Cyg OB2 #12 has been questioned in the past because 
of its advanced evolutionary state compared to the early 
O-t ype stars i n the association and its high luminosity 
llWalborn 1973T) . with the suggestion that the star might 
be a fore ground, thou gh it could also be an luminous blue 
variab le 4ciajk^_et_aU_20JL2l or a similar-luminosity binary 
(iCaballero-Nieves et al. 2014I . recently identified a compan¬ 
ion to the star, but note that it is too faint to affect the 
calculated properties of the star). 


4.5 Physical quantities derived from the HR 
Diagram 

Stellar ages and initial masses were calculated based on the 
position of each star in the H R diagram relative to the 
non-rotating Z = 0.02 models of lLeieune fe Schaerer f200ll . 
commonly used for previous studies of Cyg OB2) and the 
rotating Z = 0 .014 (the revi sed solar metallici ty value, 
lAsplund et al. 200511 models of lEkstrom et al. f2012h . both 
from the Geneva group. The rotating models have an initial 
rotation rate of Vini/v cr it = 0.4, which corresponds to typ¬ 
ical m ain s eq uence velociti es of 110-220 kms -1 for massive 
stars (Ekstrom_e£_ak_^012:), in agreement with observations 
fe.g.. [Pufton et al. 20061 7 However it is not currently clear 
whether stars are born with a narrow or broad range of ini¬ 
tial rotation velocities. 

Stellar mass loss is an important ingredi¬ 
ent in stella r models, highly influenced by the 

metallicity _ IjMe ynet et al. 19941) , rotation rate 

dMaeder fc Mevnet 2000lj. and the e xact mass-loss pre¬ 
scription used (e.g., Vink_et_aL_2001 ). The non-ro tating 
evolutionary models of Leieune &; Schaerer ('20 oHl used 
the radiative mass-loss rates from Ide Jager et al. (19881 ) 
for stars in the blue part of the HR diagram, while 
lEkstrom et al. 120121') u pdate this with the prescription 
from I Vink et al. 120011 1 where possible, and also take 
into account supra-Eddington and mechanical mass-loss 
processes in their models. Different mass-loss rates can play 
a significant role in the evolution of massive stars, their 
positions in the HR diagram, and therefore the physical 
parameters derived in such a way. These differences are 
most pronounced fo r red supergiants and the most massive 
stars (40-120 Mg, lEkstrom et al. 2012h . with the latter 
having considerably increased main-s equence lifetime s 
when mass-loss rates are increased (IMevnet et al. 1994| j. 
For moderately massive stars (20-40 M n) the effec ts of 
mass-loss can be less notable (e.s.. IMevnet et al. 19941 . note 
that the main-sequence lifetime of 25 M@ stars increased 
by only 5% when the mass-loss rate was doubled). 

The ages and masses derived from these models are 
listed in Table IB 1 1 with uncertainties. For the vast major¬ 
ity of objects this is a simple process and we do this on 
the assumption that all the stars are on their first journey 


across the HR diagram from hotter to cooler temperatures 
(generally a valid assumption since stars move back to hotter 
temperatures much faster and therefore it is rarer to observe 
stars during this transition). For stars that fall to the left of 
the zero-age main-sequence in the HR diagram we assume 
that their positions are due to uncertainties in their lumi¬ 
nosity or effective temperature and that these are zero-age 
main-sequence objects. 

To derive uncertainties on the age and initial mass 
of stars we run a MC simulation that takes into account 
both the observational uncertainties on Lboi and T e ff as 
well as uncertainties due to unresolved binary companions 
and an intrinsic spread in both luminosity and tempera¬ 
ture for a star of a given initial mass and age. The lat¬ 
ter simulates the spread in the HR diagram at a given 
initial mass and age caused by variations in fundamen¬ 
tal parameters such as rotation and metallicity, and was 
estimated to be 0.12 dex in log L/Lq and 0.02 dex in 
log T e ff from inspection of the spread in the upper main- 
sequence in rotating stella r models , and an extrapolation 
of the spread measured by iHouk et al. ( 1997l l in Hipparcos 
data. The effects of unresolved secondaries were modelled in 
our MC sim ulatio n using t he binary characteristics found by 
[Kiminki fe Kobulnickv ( 20121 ') for Cyg OB2: a binary frac¬ 
tion of 90% for O-type stars (we assume a slightly lower 
binary fraction of 80% for the B-type stars) and a flat mass 
ratio distribution over the range q = 0.005-1.0. The added 
luminosity uncertainty due to unresolved binary companions 
is generally small, typically adding an uncertainty of 0.1 dex, 
and therefore we do not consider more complex multiple sys¬ 
tems. The uncertainties shown in Figures [5] and [6] (for both 
individual stars and representative values for each spectral 
subtype) are the uncertainties on the calculated values of 
Lboi and T e // only and do not take into account these addi¬ 
tional uncertainties used when calculating the stellar mass 
and age. The median fractional uncertainty on the stellar 
mass is found to be om/M ~ 0.14 and for stellar age is 
cj t /t ~ 0.3 (for stars in the range 0-10 Myr). 


5 THE AGE OF CYG OB2 

Using the stellar ages derived from the HR diagram we can 
assess the age distribution of stars in Cyg OB2 and hence the 
star formation history of the region. This is shown in Fig¬ 
ure [7] using both the non-rotating and rotating evolutionary 
models to highlight the different age distribution that they 
each imply. We only show the 47 stars more massive than 
20 Mg (down to a main-sequence spectral type of ~B0.5V) 
based on our estimated completeness level and because the 
positions of stars less massive than this in the HR diagram 
do not provide an accurate indication of their ages (due to 
the closeness of the isochrones and the large uncertainty on 
T e ff for these stars). This excludes all the apparently over- 
luminous B-type stars discussed earlier. The uncertainties 
shown in Figure [7] were calculated from a Monte Carlo sim¬ 
ulation of the data analysis and binning process. 
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Age / Myrs Age / Myrs 

Figure 7. Histograms showing the distribution of derived ages for stars with reliably determined ages from the non-rotating (left) and 
rotating (right) models. The ages of all stars with M > 20 Mg are shown with a black line. In the left panel we also show in red the age 
distribution of stars with 20 st M/Mg C 40, which should be complete and unbiased in the age range of 0-4.5 Myr. In the right panel 
we show in red the age distribution of stars with 20 ^ M/Mq <1 35 that are unbiased in the age range of 0-6 Myr. la uncertainties are 
shown for all histograms derived from a Monte Carlo simulation of the data analysis and binning process. 


5.1 Comparison of age distributions from rotating 
and non-rotating models 

The age distribution of massive stars derived from the non¬ 
rotating models shows a strong peak at 2-3 Myrs with the 
majority of stars in the age range of 0-5 Myrs, in good agree¬ 
ment with previ ous studies of the massive star populatio n 
of Cyg OB2 (e.g.. [Hanson 20031 : [Comeron fe Pasauali 20121 ). 
Howev er, since a 2 0 Mg star ha s a lifetime of only ~8 Myr 
(in the lLeieune fe Schaerer 200ll . models) then our sample of 
massive stars with masses >20 Mg is only sensitive to stars 
<8 Myrs old. We therefore cannot reliably infer the star for¬ 
mation history greater than 8 Myrs ago when using these 
models. Furthermore this distribution is biased by the fact 
that the most massive stars have total lifetimes <4 Myrs. 
For example a 60 Mq star is estimated to have a total lifes¬ 
pan of ~3.5 Myrs, and therefore any such stars born more 
than 4 Myrs ago will not be seen now. To rectify this well- 
known bias we also show in Figure0the distribution of ages 
of stars in the mass range 20 ^ M/Mq ^ 40, which should 
be unbiased in the age range of 0-4.5 Myr. This presents a 
slightly different picture with a smaller peak suggesting that 
this peak age is an observational bias caused by the short 
lifetimes of the most massive members of the association. 

The rotating stellar models of lEkstrom et al. (20121 ) 
(Figure 0 show a different picture with a prominent but 
wide peak at an age of 4-5 Myr. The increase in the typical 
age of massive stars in Cyg OB2 when using the rotating 
stellar models can be attributed to the changes in evolu¬ 
tionary timescales induced by ro tation, which increas es the 
main sequence lifetime by ~25% dEkstrom et al. 20121 ). The 
marked influence of this on stellar ages is evident in the two 
HR diagrams (Figures [5] and [6]). For example, the group of 
04-7 supergiants with L/Lq ~ 5 x 10 s Lq have ages of 1- 


3 Myr according to the non-rotating models, but have ages 
of 3-5 Myr from the rotating models. The older ages are in 
better agreement with those derived from lower-mass stars 
in Cyg OB2, e.g., iDrew et al. f20081) estimate an ag e of 5— 
7 Myrs from A-type stars, while lWright et al. (2010bl) derive 
an age of 3.5-5.25 Myrs from an X-ray selected sample of 
approximately solar-mass stars. 


The increased stellar lifetimes of rotating stellar models 
means that the mass range over which our sample is unbi¬ 
ased has also changed relative to the non-rotating models. 
Figure 0 shows the age distribution of stars in the mass 
range 20 ^ M/Mq r 35, which should be complete in the 
age range 0-6 Myrs. This distribution is remarkably flat over 
the unbiased age range, suggesting a relatively constant level 
of star formation over the last 6 Myrs. The peak stellar age 
previously observed at 4-5 Myrs is mostly due to stars in the 
HR diagram (Figure 0) at masses of 40-60 Mq (the 04-7 
supergiants noted earlier) that have maximum stellar ages 
of <5 Myrs. The large number of slightly lower mass (25- 
50 Mq) stars with ages of 5-6 Myrs (the 07-9 supergiants 
observed in the HR diagram) could suggest that many of the 
more massive members have already evolved to their end 
states and their loss biases the observed age distribution. 


The comparisons with both sets of evolutionary mod¬ 
els indicate a number of stars in Cyg OB2 with very young 
ages, ^1 Myr. The existence of such young stars is not im¬ 
possible but appears unli kely given that long wavelength 
images (e.g., Figure 0 and lSchneider et al. 20061) show that 
the majority of dust and molecular gas been evacuated from 
the cavity in which Cyg O B2 resides. Furtherm ore both 
IWright et al. (2010bl ) and iGuarcello et al. (20131) found a 
paucity of Class I or heavily embedded low-mass stars within 
the association, suggesting the low-mass stellar population 
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is not so young (though this could also be due to the 
high levels of feed back in the vicinity o f the association, 
IWright et al. 201^ : iGuarcello et al. 20141 '). The typical un¬ 
certainties in the HR diagram shown in Figures [5] and [6] 
suggests that the positions of these stars are not incompati¬ 
ble with them being slightly older (>l-2 Myr) and the lack 
of very massive stars (>30 M@) with ages <1 Myr (where 
the large separation of isochrones allows young ages to be 
estimated more accurately) supports this idea. 

To conclude, our unbiased age distribution derived from 
rotating stellar models suggests a relatively constant level of 
star formation over the past ~6 Myr (Fig [7]) with some evi¬ 
dence for a peak in the star formation intensity between 4- 
5 Myr ago. We cannot be certain from this sample when star 
formation in the region began, but there is some evidence 
that it started at least ~6-8 Myr ago. Addressing the ques¬ 
tion of when star formation ceased is harder, but our best 
estimate is that it ended ~1 Myr ago. As will be discussed in 
Section [0. 1 1 these results are in approximate agreement with 
our star formation history fits to the mass function of stars in 
Cyg OB2. The use of different mass-loss prescriptions is not 
expected to significantly affect the evo lution of stars in ou r 
‘unbiased’ mass range of 20-35 Mq llMevnet et al. 1994ll . 
and therefore this is unlikely to explain the age spread ob¬ 
served. It is possible that the large spread in the HR dia¬ 
gram could instead be due to a single-aged population with 
an extremely large range of rotation velocities, which has 
been shown to lead to large spreads in the HR diagram 
dBrott et al. 201 ll 'l and could reproduce the entire spread 
seen in our HR diagram. However, for this to be the case 
the stars nearest the zero aged main-sequence in the HR di¬ 
agram would have to be rotating very rapidly with speeds of 
400-450 km/s, inclu ding Cyg OB2 #20 and #29, for which 
iHerrero et al. ( 19991 ') measured V sini = 25 and 180 km/s, 
respectively. This would therefore require Cyg OB2 #20 to 
have been observed virtually pole-on, which makes this sce¬ 
nario unlikely, though further high resolution spectroscopy 
of other stars in Cyg OB2 would be useful to test this hy¬ 
pothesis in more detail. It is also worth noting though that 
this scenario is not supported by measurem ents of rotation 
rates in other clusters (e.g.. iDufton et al. 20061 1. 

The uncertainties in the stellar ages we have calculated 
could explain some of this spread, but do not appear to 
be sufficiently large enough to explain the entire ~5 Myr 
spread. Furthermore the age distribution of our ‘unbiased’ 
sample shown in Figure [7] does not resemble a single-aged 
population broadened by observational uncertainties (which 
should have a more centrally-concentrated peak) and does 
appear to represent a spread of ages amongst the OB star 
population. 

5.2 Implications for large age spreads in OB 
associations 

Large age spreads have been reported in many star clus¬ 
ters and associations, typically diagnosed from luminosity 
spreads in colour-magnitude diagrams for low-mass pre- 
main sequence stars. Some authors have suggested that 
such spreads in luminosity do not equate to age spreads 
and may instead be due to intrinsic. lumi nosity spreads 
amongst a coeva l population (e.g., iHillenbrand et al. 2008l : 
iJeffries et al. 201111 potentially caused by different accretion 


histories (e.g., iBaraffe et al. 2 0091. Age spreads diagnosed 
from post-main sequence massive stars are not believed to 
suffer from such issues since the intrinsic spread and mea¬ 
surement uncertainties in Lbol and T e // are typically smaller 
than the separation of isochrones (see also IMassev 20031 1. 
Therefore if the HR diagram is sufficiently well sampled to 
overcome statistical and observational uncertainties (as it is 
here), age spreads derived from high-mass post-MS stellar 
evolution models can provide more accurate age diagnostics 
than low- and intermediate-mass pre-MS stellar evolution 
models. 

The reported age spreads for low-mass stars are often 
argued as representing a continuous period of star forma¬ 
tion within a single c luster or associatio n and in a small 
area of space (e.g., IPalla fc Stahler 20001 1. Cyg OB2 how¬ 
ever occupies a much larger region than typical young com¬ 
pact star clusters such as the Orion Nebula Cluster, with 
the stars spread over >10 pc. It has recently been argued 
that Cyg OB2 has always been a large and low density OB 
association and was never a dense and compact star cluster 
(IWright et al. 2014bl . see also lParker et al. 20141 '). It is there¬ 
fore possible that the observed age spread is due to either an 
age gradient or a series of discrete star forming events that 
have since mixed. This would imply that Cyg OB2 was not 
born as a single star cluster but as a ‘distribution’ of smaller 
groups or clusters of stars with a range of stellar ages. We 
searched for evidence of a spatial variation in the stellar age 
distribution to support this theory, but could not find any 
such evidence. If the association is composed of multiple 
populations of different ages then they are well mixed. 


6 THE UPPER END OF THE IMF IN 
CYG OB2 

In this section we analyse the IMF of the massive stars 
in Cyg OB 2, as derived usin g rotating stellar models (be¬ 
cause IWeidner fe Vink 2010 . find that such masses are in 
much better agreement with dynamical masses than those 
derived from non-rotating stellar models). Figure [S] shows a 
histogram of the masses of all 167 stars in our sample, show¬ 
ing a rapid increase as one goes from high to low masses, 
with the data flattening off for stars less massive than our 
estimated completeness limit of ~20 Mq. To calculate the 
power-law slope of the mass function we fitted this data us¬ 
ing an ordinar y least sq uares regression of Y on X (OLS 
1Y|AY llsobe et al. 1990T) of the form dN/d\ogm oc m ~ F , 
and for the mass bins >20 Mq. Because the uncertainty on 
this fit is just a fitting uncertainty we also performed a MC 
simulation to determine the true uncertainty, varying the 
masses of the individual stars according to their uncertain¬ 
ties and adjusting the number of stars in each bin as appro¬ 
priate. The la spread in the best fitting slope that results is 
then combined with the fitting error to produce the true un¬ 
certainty on the mass function slope. The best fitting mass 
function slope for all stars in our sample is V = 2.06 ± 0.06, 
significant ly steeper than the canonical Salpeter slope of 
T ~ 1.35 (lSalpeter 1955! ). although it agrees with that de- 
termi ned by other authors for the massive stars in Cy g OB2 
fe.s.. iKiminki et al. 20071 : IComeron fe Pasquali 20121 . Note 
however that this uncertainty is still falsely low because it 
does not take into account the random sampling of the IMF. 
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Figure 8. Mass functions for massive stars in Cyg OB2 derived using the rotating stellar models of lEkstrom et al. 120121 . Left: Histograms 
showing the distribution of derived masses for all 167 stars in our sample (black histogram, with ltr error bars determined from MC 
simulations) and for an age-unbiased sample of 82 stars aged <3.5 Myrs ( red histogram, wi th 1 a error bars). Both histograms have been 
fitted using an ordinary least squares regression of Y on X (OLS (FIXI. Ilsobe et al. 199011 for the mass bins >20 Mg, giving slopes of 
r = 2.06 ± 0.06 and 1 = 1.39 ± 0.19 for the complete and age-unbiased samples respectively (note that these uncertainties do not take 
into account random sampling of the IMF). Right: A cumulative mass function showing N (M* < M) as a function of M for all 98 stars 
in our sample more massive than 10 Mg (black line) with the la uncertainty (calculated from MC simulations) shown in grey. Also 
shown are the results of simulated IM Fs for T = 1.3 (red line, with la random sampling uncertainty shown in pink) and F = 2.2 (green 
line) calculated using the equations of lMaschberger F201 4 see full explanation in the text). The blue line shows a simulated and evolved 
IMF for T = 1.3 for stars aged 3—7 Myr (i.e. a constant star formation rate from 3 to 7 Myr ago, and no star formation since then). The 
uncertainties for three simulated IMFs are similar and so only one is shown for clarity. 


The mass distribution of stars in our sample is likely to 
be altered by the loss of the most massive stars that have al¬ 
ready evolved to their end states (Section O, steepening the 
observed mass function. To counter this we have attempted 
to identify an unbiased sample of stars from which to study 
the initial mass function in Cyg OB2, selecting all stars with 
ages <3.5 Myrs, which we estimate from the rotating stel¬ 
lar models (Figure [6j) as the age at which the most mas¬ 
sive stars in Cyg OB2 will have evolved to their end states. 
The mass function of this unbiased sample of 82 stars is 
shown in Figure [8l which shows a slightly shallower slope 
to that of the full sample. An OLS (V|-Y) fit gives a slope 
of T - 1.39 ± 0.19, in good agreement with the concept of 
a universal IMF at high masses with a slope of T = 1.3 
llKroupa et al. 200ll ; iBastian et al. 20ld ). This value also 
agrees with the mass function measu red at lower masses 
(r = 1.09 ± 0.13, IWright et al. 2010bl), and the mass func - 
tion of massive stars found by Massey fe Thompson f!99ll ) 
of r = 1.0 ± 0.3, although their sample was considerably 
smaller than the current sample and may therefore have 
been incomplete. 


6.1 The Cumulative Mass Function 

To avoid the well-known and negative effects of binning our 
data, and to overcome any uncertainties arising from the 
narrow age range used for the unbiased mass sample we 
have also considered the mass function of stars in Cyg OB2 
represented as a cumulative function (CF). This is shown in 
Figure [8] with 1 <t uncertainties calculated from a MC simu¬ 


lation using the full uncertainties on all stellar masses (the 
la uncertainties shown are very small due to the cumulative 
nature of the function and the lack of binning). 

The CF shows a steady increase up to ~30 Mg, and 
then flattens off, reaching a total of 98 stars with masses 
> 10 M 0 . We compare our CF with si mulated IMFs gener¬ 
ated using the equations described in IMaschberger (20121) . 
varying a (= F + 1), but fixing /3 = 1.4 (the low-mass expo¬ 
nent has no effect on the high-mass IMF, but is required for 
the calculations) and drawing individual stellar masses from 
the mass range 0.01-150 Mg. We continue to do this until 
the number of stars with M > 20 Mg reaches the observed 
number of 98 and use this to produce a cumulative function 
scaled to the observed number of stars with M < 20 Mg. 
The simulated IMF is generated 10,000 times and the me¬ 
dian CF shown in Figure [5] with ±1 <t range illustrated. 

We initially simula ted an unevolved ( i.e. counting all 
stars) ‘universal’ IMF (|Kroupa et al. 20011 ) with T = 1.3, 
but found that this over-predicted the number of very mas¬ 
sive stars (> 50 Mg) relative to the moderately massive stars 
(20-50 Mg) and we can reject a standard, unevolved mass 
function of this type with a confidence of ~3<r. Varying T 
we found a best fit with an un-evolved IMF for T = 2.2+0 ®, 
in approximate agreement with that found fitting a binned 
IMF for our entire sample (but with larger uncertainties that 
also take into account the random sampling of the IMF). 

We then simulated a number of evolved mass functions 
using the ‘universal’ T = 1.3 IMF and different star forma¬ 
tion histories. For a star of a given mass and age we com¬ 
pute its position in the HR diagram according to the models 
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of lEkstrom et al. 120121 ) and remove those stars that have 
evolved beyond their reddest position in the HR diagram. 
This excludes WR stars from the simulated mass functions, 
matching our observed CF that also leaves them out (be¬ 
cause of the difficulty of determining their initial masses). 

We found that a number of different star formation his¬ 
tories provided satisfactory fits (within lcr) to the observed 
mass function, preventing us from accurately determining it 
in this way, though we were able to place some partial con¬ 
straints on it. We find fits using a ‘universal’ T = 1.3 IMF 
requires the star formation rate over the last ~3 Myr to be 
very low, <20% of the total star formation in Cyg OB2, with 
the fit improving slightly as the fraction of star formation 
in this time period decreases. Exploring models with this 
feature we find acceptable fits when the duration of star 
formation is at least 2 Myrs (i.e., a constant level of star 
formation from 3 to 5 Myrs ago), with the best fit obtained 
for a model with a constant star formation rate over the 
time period of 3 to 7 Myr ago, in reasonable agreement with 
that found from the distribution of stellar ages studied ear¬ 
lier. Significant star formation prior to ~7 Myr ago degrades 
the fit to the cumulative mass function. These results are in 
good agreement with the stellar age distribution discussed 
in in Section [5] 

To conclude, the current mass function of massive stars 
in Cyg OB2 is relatively steep, with a slope of F ~ 2.1, 
but the stellar age distribution suggests that it has been 
steepened by the most massive stars evolving to their end 
states. From both an unbiased sample of stars and a mod¬ 
elling an evolved cumulative mass function we find that the 
mass functio n is consistent with the canonical T =1.3 ‘uni¬ 
versal’ IMF (ISalpeter 19551 : iKroupa et al. 20oH ). 


6.2 Have their been any supernovae in Cyg OB2? 

There is considerable evidence, both in the age distribution 
of stars and in the steepened IMF, to suggest that some of 
the most massive stars in Cyg OB2 have already evolved to 
their end states and therefore that there has already been a 
supernova within the region. 

Despite this there are not any know n supernova r em- 
nants (SNRs) within Cyg OB2, although iGreen f 2009h do 
list a number of SNRs in the wider Cygnus X region 
(including 7 -Cygni, G78.2+2.1) but these are predomi¬ 
nantly to t h e we st in the Cyg OBI and OB9 associations. 
iButt et al. ('200311 identified an expanding shell of molecu¬ 
lar and dusty material from CO J = 1 —)• 0 and IRAS 
observations that they suggested could be due to a SNR, 
but this cavity could equally arise from the combined ef¬ 
fects of ionising radiation and stellar winds from the massive 
stars in Cyg OB2. The lack of observed SNRs in Cyg OB2 
does not however rule out any previous supernovae in the 
association, since SNRs t ypic ally remain visible for only 
~10,000 yrs le.g.. ILevenson et al. 199a) and may not leave 
any detectable signature when expanding into hot, low- 
density me dia suc h as the cleared-out cavities in OB as¬ 
sociations (IChu 19971 '). 

There have been a number of very high-energy sources 
identified towards Cyg OB2 that have been cons i dered 
as possible comp act stellar remnants (e.g., iBednarek 20031 : 
IButt et al. 2006I ). but confirming their nature or identi¬ 
fying counterparts at other wavelengths has proved diffi¬ 


cult. |Abdq^et^al~(200^) discovered a ~115,800 yr old 7 -ray 
pulsar (J2032+4127) with a frequency of 6.98 Hz toward 
Cyg OB2, which ha s since been detect ed in the radio with 
the same frequency dCamilo et al. 20091) and is also believed 
to be associated with th e previousl y-unidentified TeV source 
J20 32+4130 (lAh aronian et al. 20021) as its pulsar wind neb¬ 
ula. ICamilo et al. 120091 ) argue that the pulsar is likely to 
be located within the Cyg OB2 association, and combined 
with our newly derived age of ~l-7 Myrs this is consistent 
with the pulsar originating within Cyg OB2 supporting the 
theory that the association has already seen it’s first super¬ 
nova. 

Another indication of past supernovae is the presence of 
high velocity ‘runaway’ OB stars that are eje cted when their 
binary com panion expl odes as a supernova dBlaauw 19611) . 
IComeron fc Pasquali (20071) identified an 04If runaway star 
with a proper motion of ~40 km/s moving away from 
Cyg OB2 and argued that the star was likely ejected from 
the association ~1.7 Myrs ago. The authors estimate a mass 
of 70 ± 15 M@ and an age of ~1.6 Myrs (by comparisons 
with non-rotating stellar isochrones) suggesting the star was 
ejected from Cyg OB2 shortly after it was born, and there¬ 
fore not from a supernova explosion. The other mechanism 
for stellar eje ction is via dynamical encounters in dense stel¬ 
lar systems (|Poveda et al. 19671) . yet this process requires 
much higher stellar densities t ha n are believed to have ex¬ 
isted in Cyg OB 2 (IWright et al. 2014bl) making it unlikely . 
Using the rotating stellar models of lEkstrom et al. f20121) . 
and with updated photometry we re-calculated the proper¬ 
ties of this star. We find B — K s = 5.27 mag (using Tycho- 2 
and 2MASS photometry. iHpg et al. 2000 l : ICutri et al. 20031 ). 
which for a star of type 04If with ( B — K s )o = —1.15 mag 
dMartins fc Plez 20061 ) gives Ay = 5.22 mag (using the red¬ 
dening law determined above) and therefore log Lbol = 
5.85 L 0 and T e // = 40,422 K. This implies an initial mass 
of ~55 Mq and an age of 4-5 Myr. This would suggest the 
star was ejected from Cyg OB2 when it was about ~3 Myr 
old, the age at which the most massive star in a population 
might explode as a supernova. This is therefore consistent 
with a scenario whereby the star’s very massive (>100 M 0 ) 
companion exploded as a supernova and ejected it’s com¬ 
panion from the association. The remnant of the massive 
star was presumably ejected from Cyg OB2 in the opposite 
direction and could be searched for to verify this scenario. 


7 THE TOTAL MASS OF CYG OB2 

Using a sample of massive stars that is believed to be com¬ 
plete within a given mass range, and the assumption that 
the IMF has a ‘universal’ form in all star clusters and OB 
associations, it is possible to estimate the total mass of a 
group of stars beyond that which have been detected. To 
do this for Cyg OB2 we have selected a mass range of 20- 
40 M 0 for which we believe our sample is complete (at least 
for stars aged <6 Myrs). Above 40 M 0 our sample will be 
incomplete for stars aged <6 Myrs and below 20 M 0 we 
have argued that our sample is likely to be observationally 
incomplete. Based on the stellar masses determined above 
the number of stars in this mass range is 36^4 (uncertainties 
calculated from MC simulations). 

To estimate the total mass of Cyg OB2 we ran a MC 
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simulation by sampling stellar masses fr om a ‘universal’ IMF 
with a = 2.30 (r = 1.3) and /3 = 1.40 llKroupa et al. 200 ll : 
iMaschberger 20121 4 across the entire mass range. We pro¬ 
ceeded to count the total mass of stars produced and the 
number of stars in the mass range of 20-40 M 0 , until the 
latter number reached the observed number of stars in that 
mass range. We performed this simulation 10,000 times, 
recording the total stellar mass that was necessary in each 
iteration to produce the observed number of stars. Further¬ 
more we varied the observed number of stars, 36 I 4 , in each 
iteration according to the uncertainties on that number. 
This MC simulation therefore takes into account both the 
uncertainty in the total number of stars observed, and the 
variation in the number of stars in a cluster of a given mass 
due to the random sampling of the mass function. 

We find that the observed number of massive stars in 
Cyg OB2 in the mass range of 20-40 M@ can be reproduced 
if the association has a total stellar mass of 16500l2goo M 0 . 
This is based on three assumptions, firstly that we are com¬ 
plete in the chosen mass range, as argued above. Secondly, 
that the vast majority of star formation in Cyg OB2 oc¬ 
curred in the last 6 Myr (broadly supported by the age 
distribution and the mass function simulations, though it 
is difficult to test this assumption without extending the 
completeness of the sample to lower masses). The final as¬ 
sumption is that the IM F in Cyg O B2 is we ll represented 
by the ‘universal’ IMF of lKroupa et al. (200111 . This is sup¬ 
ported by our IMF simulati ons at high stellar masses and 
that found at lower masses llWright et al. 2010 bT ). however 
a small variation can lead to large differences in the total 
stellar mass. For example increasing (decreasing) the high- 
mass exponent a by 0.1 increases (decreases) the total stel¬ 
lar mass by 5600 (3900) M 0 . However, for comparison with 
other young clusters and associations that have also adopted 
such an IMF we believe this is the most appropriate choice. 

Our mass estimate is smaller than s ome previous mass 
estimates of Cyg OB2, notably that of IKnodlseder f2000l i 
who esti mated a mass o f (4 — 10) x 10 4 M 0 for the as¬ 
sociation. IKnodlseder f 200C1 I based their estimate on the 
number of stars brighter than a given magnitude (equiva¬ 
lent to an F3V star) from a background-subtracted 2MASS 
star counts study. The choice of a background region in this 
crowded and complex area of the Galactic Plane is fraught 
with difficulties and could lead to large uncertainties in any 
resulting star counts or mass estimates. IKnodlseder ('20001 1 
also estimate the total number of O stars in Cyg OB2 to 
be 120 ± 20 , an estimate almost double the size of our cen- 
sus and one which other studies have not supported (e.g., 
lHanson 2003lh Our mass estimate is more consistent with 
recent estimates bas ed on studies of lower-mass stars, e.g., 
IWright et al. f2010bll estimate d (2 —4) x 10 4 M 0 from an X- 
ray study and Drew et al. 1200^) estimated (1 —4) x 10 4 M 0 
from a photometric study of A-type stars in the association. 


8 CONCLUSIONS 

We have compiled a census of all known massive stars in the 
nearby OB association Cygnus OB2. We have gathered data 
from across the literature to achieve a census containing 169 
primary OB stars, including 52 O-type stars and 3 WR stars. 
We have used the sample to measure the extinction 


l aw towards Cyg O B2, fittin g the available photometry to 
a lFitzpatrick feMassaf 20071 1 Rv = 2.91 ± 0.06 extinction 
law. This is in good agreement with that determined by 
lHanson 120031 ) in the optical and differing only in the near- 
IR due to shifts in the intrinsic near-IR colours of O-type 
stars in recent years. We use this new law to derive individ¬ 
ual extinctions for all the stars in our sample, the majority 
of which have extinctions Ay = 4-7 mag, with a typical 
uncertainty of ~0.3 mag. 

A Hertzprung-Russell diagram was compiled and com¬ 
pared to both non-rotating and rotating stellar evolution 
tracks, from which stellar masses and ages are calculated. 
Uncertainties are derived for all quantities using a MC 
simulation that takes into account all available uncertain¬ 
ties including quantified classification uncertainties, photo¬ 
metric uncertainties, extinction law uncertainties, estimated 
spreads in fundamental parameters, and unresolved binaries. 

The distribution of stellar masses is assessed using both 
a binned mass function and a cumulative mass function, 
both of which are affected by the loss of the most massive 
stars to stellar evolution. The current mass function is rela¬ 
tively ste ep with a slope of T = 2.2 (in agreement with that 
found bv lKiminki et al. 2007f) . though when considering ei¬ 
ther an unbiased sample of stellar masses or by modelling 
the effects of stellar evolution on the mass function we find 
it to be consistent with a Salpeter-like ‘universal’ slope of 
T = 1.3. 

The non-rotating stellar evolution models suggest a 
peak stellar age of 2-3 Myr, as derived by previous authors, 
while the rotating stellar models suggest a slightly older 
peak age of 4-5 Myr. The results from the rotating stellar 
models are favoured, both for physical reasons, and because 
the age distribution that they imply is in better agreement 
with that found from studies of low- and intermediate-mass 
stars. When considering only stars in a mass range unaltered 
by evolutionary effects the age distribution is not inconsis¬ 
tent with a relatively constant level star formation over at 
least the past ~6 Myrs. Our star formation history fits to 
the mass function of stars in Cyg OB2 suggest star forma¬ 
tion most likely began ~7 Myr ago. The lack of molecular 
material or dust in the immediate vicinity of the associa¬ 
tions suggests star formation likely ceased at least Myr 
ago, and this is also supported by our mass function mod¬ 
els. Our best estimate is that the majority of star formation 
occurred during the period of 1-7 Myr ago, with a possible 
in the star formation intensity between 4-5 Myr ago. 

Given the evidence from both the age distribution and 
the mass function that some of the most massive stars have 
already evolved to their end states we searched the litera¬ 
ture for evidence of this and found both a young pulsar and 
a runaway O star that may have originated in Cyg OB2, 
consistent with the view that Cyg OB2 has already seen its 
first supernova. 

Finally we use a MC simulation to estimate the to¬ 
tal mass of the association, taking into account the loss of 
the most massive stars, and estimate a total ste llar mass of 
16500 ~t 28 oo M 0 (assuming a lKroupa et al. 200 ll . IMF) that 
has formed over the past 6 Myr. This is at the lower end of 
many previous estimates, based on the number of low- and 
intermediate-mass stars, but consistent given the uncertain¬ 
ties involved in extrapolating the IMF. 

I 11 conclusion our new census of massive stars in 
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Cyg OB2 provides an updated view of the high-mass stel¬ 
lar population in one of the largest groups of young stars in 
our Galaxy and for the first time these results show a good 
agreement with those determined for low- and intermediate- 
mass stars, thereby providing a coherent picture for this 
massive OB association. In a future paper we will use new 
photometry and spectroscopy to search for previously undis¬ 
covered massive stars, extending this census to both lower- 
luminosity B-type stars in Cyg OB2, and O and B-type stars 
over a wider area and further along the sightline. 
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APPENDIX A: UVEX PHOTOMETRY AND 
TRANSFORMATIONS TO THE UBV SYSTEM 

iMassev & Thompson (T99lh note that a number of their 
photometric measurements suffer from blended photometry, 
of which 11 are included in our census. To determine the 
severity of this blending and to provide alternative photom¬ 
etry for these sources we gathered Ug' r' photometry from 
UVEX, the Ultra-Violet Excess Survey (iGroot et al. 2009h . 
selecting sources in the magnitude range 13-19 in each fil¬ 
ter. The Ug'r' filters are directly transferable onto the UBV 
photometric system using transformations presented on the 
CASU websitc0. However these transformations were de¬ 
rived from observations of Landolt standard fields, which 
are typically dominated by unreddened late-type field stars 
that can have different spectral energy distributions to 
the reddened early-type stars in our sample. We there- 
fore calculated our own tra nsformations by comparing the 
iMassev fe Thompson (T99li l UBV photometry with UVEX 
Ug'r' photometry and fitting new transformations (exclud¬ 
ing the blended stars). Figure [All shows the results of the 
least-squares straight-line fits from which we derive the 
transformations 

UuVEX — U Lan = —0.22 — 0.180(1/ — B)Lan (Al) 

g — B = 0.024 — 0.523(£ - V)Lan (A2) 

2 Cambridge Astronomical Surveys Unit: Colours and Trans¬ 
formations for the Isaac Newton Telescope Wide Field Survey: 
http: //www.ast.cam.ac.uk/~wfcsur/technical/photom/colours/ 
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Figure Al. Calibration of transformations between thelMassev &; Thompson f 199 ij) UBV photometry and the UVEX Ug’r’ photometry. 


All stars are shown as black dots, with the stars listed bv lMassev fe Thompson i 1 99 if) as ‘blended’ shown with red dots. The best fitting 
transformation obtained using a least squares fit that excludes the ‘blended’ stars are shown with red dashed lines and are listed in each 
panel with the number of stars used for the fit and the dispersion on the best fit. 


(U — g')uvEx = 0.544 + 1.303(17 — B)Lan (A3) 

g -r = 0.213 + 0.855(B - V) Lan (A4) 

which are based on 101, 96, 90 and 60 stars respectively, and 
have dispersions of 0.080, 0.068, 0.056, and 0.058 mag re¬ 
spectively (no stars were found further than 3<r from the best 
fits). These are in good agreement with the existing CASU 
transformations, differing from them by less than 10% in the 
slope and typically by ~0.1 in the intercepts. 


APPENDIX B: THE FINAL CATALOG 

The full census of known OB stars in Cyg OB2 is presented 
in Table [BT] in order of increasing primary spectral sub- 
type and with derived observational (Av), physical (Te//, 
Lto(), and stellar (M, age) qua ntities (the latter calculated 
using the lEkstrom et al. 20121 . rotating evolutionary mod¬ 
els). Stars that fall to the left of the zero-aged main se¬ 
quence in the HR diagram are listed as having ages of 0 Myr, 


while those whose positions in the HR diagram suggest ages 
>10 Myr are listed as such (see Section [4.51 for discussion 
of these issues). References for spectral classifications are 
noted at the end of the table. The full table, which also con¬ 
tains UB V and J HK S p hotometry, st ellar quantities derived 
using the iLeieune fe Schaerer f 200lh non-rotating models, 
along with 1 <t uncertainties on all derived properties, is avail¬ 
able from Vizier. 


© 0000 RAS, MNRAS 000, 000-000 














20 Wright et al. 


Table Bl: Massive stars in Cygnus OB2 


RA 

(J2000.0) 

Dec 

(J2000.0) 

Spectral 

type 

Spectral 

reference 

S58 

Source number 

MT91 C02 

Ay 

(mag) 

log L 
( l q) 

T eff 

K 

M 

(M e ) 

Age 

(Myrs) 

20:32:03.10 

41:15:19.9 

WC4 

vH88 



WR144 

- 

4.99 

5.05 

>20 

2-8 

20:32:06.28 

40:48:29.7 

WN7o/CE+07V((f)) 

S96,M09 



WR145 

- 

5.57 

4.70 

>25 

2-7 

20:35:47.09 

41:22:44.7 

WC 6 +O 8 III 

DOO,DOO 



WR146 

7.83 

5.06 

4.90 

>25 

2-4 

20:33:08.78 

41:13:18.1 

03H+06V 

W02,W02 

22 

417 


6.92 

5.82 

4.63 

49.9 

4.12 

20:33:14.16 

41:20:21.5 

03If 

MT91 

7 

457 


5.30 

5.63 

4.63 

46.7 

2.67 

20:33:18.02 

41:18:31.0 

05111 

K07 

8C 

483 


4.85 

5.57 

4.61 

41.6 

3.34 

20:33:10.74 

41:15:08.0 

051+03.5III 

NZ12,NZ12 

9 

431 


6.76 

5.85 

4.59 

51.6 

4.26 

20:34:08.54 

41:36:59.3 

O5H+B0V 

K12a,K12a 

11 

734 


5.40 

5.73 

4.59 

43.7 

4.56 

20:33:23.46 

41:09:12.9 

05.5 V 

MT91 


516 


7.43 

5.85 

4.60 

51.6 

4.21 

20:33:13.25 

41:13:28.6 

06V 

MT91 


448 


7.22 

5.03 

4.59 

28.6 

0.76 

20:31:37.50 

41:13:21.1 

06IV+09III 

K08,K08 

3 



5.84 

5.54 

4.59 

38.0 

4.18 

20:33:15.18 

41:18:50.1 

061+05.5III 

dB04,dB04 

8A 

465 


4.84 

5.70 

4.57 

41.1 

5.00 

20:33:14.84 

41:18:41.4 

06.5111 

MT91 

8B 

462 


5.12 

5.48 

4.57 

35.2 

4.56 

20:34:44.10 

40:51:58.0 

06.5111 

N08 



A24 

7.29 

5.19 

4.57 

29.6 

3.41 

20:33:40.88 

41:30:18.5 

07V 

MT91 


611 


5.54 

4.66 

4.57 

22.3 

0 

20:34:13.50 

41:35:02.6 

07V 

MT91 

29 

745 


5.36 

4.94 

4.57 

25.4 

1.78 

20:34:29.52 

41:31:45.5 

07V+09V 

K08.K08 


771 


6.93 

5.18 

4.57 

29.0 

3.50 

20:31:59.61 

41:14:50.4 

07V 

C12 



E45 

5.89 

4.79 

4.57 

23.3 

0.17 

20:32:13.77 

41:27:12.7 

07111 

K07 

4 

217 


4.39 

5.20 

4.56 

28.6 

4.14 

20:32:22.43 

41:18:19.0 

07I+06I+09V 

B76,B76,D13 

5 



6.02 

6.20 

4.54 

93.1 

2.95 

20:31:36.91 

40:59:09.1 

071 

N08 



A15 

8.26 

5.42 

4.54 

31.8 

5.25 

20:30:27.30 

41:13:25.0 

071+091 

S10.S10 



B17 

8.08 

5.03 

4.54 

24.8 

4.17 

20:33:17.49 

41:17:09.2 

07.5V 

MT91 

24 

480 


5.63 

5.02 

4.55 

25.4 

3.40 

20:33:26.77 

41:10:59.5 

07.5V 

MT91 


534 


6.36 

4.84 

4.55 

23.4 

1.77 

20:32:31.50 

41:14:08.0 

07.5111 

K12a 


267 


8.03 

5.51 

4.56 

34.7 

5.17 

20:31:10.57 

41:31:53.0 

08V 

K08 

1 

059 


5.18 

5.08 

4.54 

25.5 

4.46 

20:32:27.67 

41:26:21.7 

08V 

MT91 

15 

258 


4.36 

4.79 

4.54 

21.9 

2.45 

20:32:38.58 

41:25:13.6 

08V 

H03 

16 

299 


4.40 

4.92 

4.54 

23.4 

3.65 

20:32:45.45 

41:25:37.3 

08V 

MT91 

6 

317 


4.54 

5.04 

4.54 

24.8 

4.28 

20:32:59.17 

41:24:25.7 

08V 

MT91 


376 


4.90 

4.69 

4.54 

20.9 

1.04 

20:33:02.94 

41:17:43.3 

08V 

MT91 


390 


6.63 

4.93 

4.54 

23.5 

3.70 

20:33:13.67 

41:13:05.7 

08V 

MT91 


455 


6.13 

4.75 

4.54 

21.4 

1.89 

20:33:18.08 

41:21:36.6 

08V 

MT91 


485 


5.32 

4.79 

4.54 

21.8 

2.40 

20:33:30.43 

41:35:57.5 

08V 

MT91 

74 

555 


6.45 

5.04 

4.54 

24.9 

4.31 

20:32:34.80 

40:56:17.0 

08V 

N08 



A38 

6.31 

4.64 

4.54 

20.3 

0.41 

20:34:21.95 

41:17:01.6 

08III+08III 

K09,K09 

73 



6.27 

4.62 

4.53 

19.6 

1.36 

20:33:02.90 

40:47:25.0 

0811 

H03 



A20 

8.13 

5.58 

4.53 

35.0 

5.72 

20:32:50.03 

41:23:44.6 

08.5V 

MT91 

17 

339 


4.92 

4.79 

4.53 

21.2 

3.67 

20:33:16.36 

41:19:01.9 

08.5V 

MT91 

8D 

473 


5.13 

4.69 

4.53 

20.2 

2.59 

20:33:21.04 

41:17:40.1 

08.5V 

MT91 


507 


5.41 

4.53 

4.53 

18.7 

0.13 

20:33:25.67 

41:33:26.6 

08.5V 

MT91 

25 

531 


5.46 

4.99 

4.53 

23.5 

4.88 

20:31:45.39 

41:18:26.8 

08.51 

MT91 


138 


6.57 

5.07 

4.51 

23.6 

6.11 

20:31:18.31 

41:21:21.7 

09V+B6V 

K12a,K12a 


070 


6.86 

4.55 

4.52 

18.4 

1.97 

20:32:16.53 

41:25:36.4 

09V 

MT91 

14 

227 


4.55 

4.66 

4.52 

19.1 

3.61 

20:33:09.58 

41:13:00.6 

09V+B9V 

K09,K09 


421 


6.77 

4.26 

4.52 

16.3 

0 

20:34:04.95 

41:05:13.2 

09V 

MT91 

41 

716 


6.12 

4.45 

4.52 

17.5 

0.26 

20:34:09.52 

41:34:13.4 

09V 

MT91 

75 

736 


5.43 

4.50 

4.52 

18.0 

1.07 

20:33:09.41 

41:12:58.2 

09V 

K07 


420 


6.54 

4.06 

4.52 

15.2 

0 

20:31:49.65 

41:28:26.8 

09111 

K09 

20 

145 


4.12 

4.42 

4.50 

16.8 

1.46 

20:33:46.15 

41:33:00.5 

091 

K07 

10 

632 


5.66 

5.68 

4.50 

37.4 

5.54 

20:33:15.74 

41:20:17.2 

09.5V 

MT91 

23 

470 


5.13 

4.42 

4.50 

16.8 

1.35 

20:33:59.57 

41:17:36.1 

O9.5V+B0V 

MT91,K12 

27 

696 


5.81 

4.48 

4.50 

17.2 

2.58 

20:30:57.70 

41:09:57.0 

09.5V 

N08 



A26 

6.69 

4.69 

4.50 

18.7 

5.31 

20:34:16.05 

41:02:19.6 

09.5V 

C12 



E54 

6.23 

4.54 

4.50 

17.6 

3.62 

20:31:46.00 

41:17:27.4 

09.51 

MT91 


140 


2.17 

4.54 

4.48 

16.6 

5.97 

20:32:13.07 

41:27:24.9 

BOV 

MT91 


213 


4.30 

4.18 

4.49 

14.5 

0 

20:32:59.61 

41:15:14.6 

BOV 

MT91 

41 

378 


6.93 

4.65 

4.49 

17.5 

6.55 

20:33:10.10 

41:13:10.1 

BOV 

MT91 

51 

425 


6.55 

4.44 

4.49 

16.0 

4.32 

20:33:10.34 

41:13:06.4 

BOV 

MT91 


426 




Unresolved 


20:33:10.50 

41:22:22.8 

B0V+B3V 

MT91,K12 


429 


5.45 

3.99 

4.49 

13.5 

0 

20:33:37.02 

41:16:11.4 

BOV 

MT91 

70 

588 


5.85 

4.69 

4.49 

17.9 

6.80 

20:33:59.32 

41:05:38.4 

BOV 

MT91 


692 


5.74 

4.14 

4.49 

14.2 

0 

20:30:51.12 

41:20:21.6 

BOV 

K07 


020 


7.20 

4.36 

4.49 

15.5 

2.65 
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Table Bl: Massive stars in Cygnus OB2 


RA 

(J2000.0) 

Dec 

(J2000.0) 

Spectral 

type 

Spectral 

reference 

S58 

Source number 
MT91 C02 

A v 

(mag) 

log L 
( l q) 

T eff 

K 

M 

(M e ) 

Age 

(Myrs) 

20:32:58.79 

41:04:29.9 

B0V+B2V 

K07,K09 


372 


6.98 

3.85 

4.49 

12.9 

0 

20:33:11.02 

41:10:31.9 

BOV 

K07 


435 


7.13 

4.20 

4.49 

14.6 

0 

20:33:17.40 

41:12:38.7 

BOV 

K07 


477 


6.49 

4.12 

4.49 

14.1 

0 

20:33:18.56 

41:24:49.3 

BOV 

K07 


490 


5.77 

3.76 

4.49 

12.5 

0 

20:33:42.38 

41:11:45.8 

BOV 

K07 


620 


6.24 

4.24 

4.49 

14.7 

0.19 

20:32:39.51 

40:52:47.5 

BOV 

C12 



E48 

6.60 

4.69 

4.49 

17.9 

6.78 

20:32:34.90 

40:52:39.0 

B0.2V 

N08 



A33 

6.80 

4.42 

4.47 

15.1 

6.62 

20:33:21.14 

41:35:52.0 

BOIII 

K07 


509 


6.66 

4.14 

4.49 

14.2 

0 

20:33:39.14 

41:19:26.1 

BOIab 

MT91 

19 

601 


5.72 

5.27 

4.49 

26.0 

6.74 

20:34:44.70 

40:51:46.0 

BOIa 

H03 



A27 

7.58 

5.63 

4.49 

35.2 

5.94 

20:34:58.70 

41:36:17.0 

BOIb+BOIII 

H03,K09 



A36 

6.05 

5.24 

4.49 

25.5 

6.78 

20:32:39.06 

41:00:07.8 

BOIa 

C12 



E47 

6.53 

5.83 

4.49 

42.3 

4.79 

20:33:38.22 

40:53:41.2 

BOIb 

C12 



E52 

5.91 

4.97 

4.49 

21.2 

7.06 

20:32:03.74 

41:25:10.9 

B0.5V 

MT91 


187 


5.31 

4.05 

4.44 

12.3 

3.74 

20:32:27.76 

41:28:51.9 

B0.5V 

MT91 

21 

259 


3.80 

4.20 

4.44 

12.9 

7.63 

20:33:39.84 

41:22:52.4 

B0.5V+B2.5V 

K12.K12 


605 


4.66 

4.10 

4.44 

12.5 

5.39 

20:34:06.10 

41:08:09.6 

B0.5V+B1.5V 

K12.K12 


720 


6.95 

4.20 

4.44 

12.9 

7.73 

20:32:39.50 

40:52:47.0 

B0.5V 

N08 



A31 

7.20 

4.63 

4.44 

15.7 

>10 

20:33:06.62 

41:21:13.3 

B0.5V 

K07 


409 


5.96 

3.91 

4.44 

11.6 

0 

20:33:18.70 

40:59:37.9 

B0.5III 

C12 



E51 

6.63 

4.54 

4.44 

14.9 

>10 

20:32:13.53 

41:27:30.0 

B1V 

MT91 


215 


4.05 

3.58 

4.39 

9.2 

0 

20:32:37.03 

41:23:05.1 

B1V 

MT91 


292 


5.31 

4.09 

4.39 

10.9 

>10 

20:32:38.87 

41:25:20.8 

B1V 

K07 


300 


4.34 

3.62 

4.39 

9.2 

0.90 

20:33:04.42 

41:17:08.9 

B1V 

K07 

54 

395 


5.96 

3.90 

4.39 

10.2 

>10 

20:33:05.22 

41:17:51.6 

B1V 

K07 


400 


5.54 

3.74 

4.39 

9.6 

6.58 

20:33:15.37 

41:29:56.6 

B1V 

MT91 


467 


5.41 

3.97 

4.39 

10.4 

>10 

20:33:23.24 

41:13:41.9 

B1V 

MT91 

66 

515 


6.77 

4.01 

4.39 

10.6 

>10 

20:31:33.38 

41:22:49.0 

B1V+B2V 

K07,K12a 


103 


6.37 

4.04 

4.39 

10.7 

>10 

20:32:14.56 

41:22:33.7 

B1V 

K07 


220 


5.38 

3.61 

4.39 

9.1 

0.61 

20:32:21.35 

41:18:35.5 

B1V 

K07 


238 


5.71 

3.45 

4.39 

8.7 

0 

20:32:50.69 

41:15:02.2 

B1V 

K07 


343 


6.55 

4.00 

4.39 

10.6 

>10 

20:32:56.66 

41:23:41.0 

B1V 

K07 


365 


4.81 

3.62 

4.39 

9.2 

1.14 

20:33:10.46 

41:20:57.6 

B1V 

K07 


428 


6.00 

3.95 

4.39 

10.4 

>10 

20:33:19.16 

41:17:44.9 

B1V 

K07 


492 


5.67 

3.50 

4.39 

8.9 

0 

20:33:23.37 

41:20:17.2 

B1V 

K07 


517 


5.24 

3.79 

4.39 

9.8 

9.10 

20:33:42.57 

41:14:56.9 

B1V 

K07 


621 


6.46 

3.79 

4.39 

9.8 

9.02 

20:34:04.43 

41:08:08.4 

B1V 

K07 


712 


6.15 

4.18 

4.39 

11.4 

>10 

20:34:24.56 

41:26:24.7 

B1V 

K07 


759 


5.85 

3.68 

4.39 

9.4 

3.52 

20:30:39.70 

41:08:48.0 

B0.7Ib 

H03 



A23 

7.38 

5.42 

4.42 

26.3 

7.04 

20:33:47.88 

41:20:41.7 

B1III 

MT91 

26 

642 


5.77 

4.84 

4.39 

15.9 

>10 

20:33:15.51 

41:27:32.9 

B1III 

K07 


469 


5.13 

3.81 

4.39 

9.8 

9.88 

20:33:46.85 

41:08:01.9 

B1III 

K07 


635 


5.78 

3.96 

4.39 

10.4 

>10 

20:31:22.03 

41:31:28.0 

B1I 

MT91 

2 

083 


3.98 

4.61 

4.39 

14.3 

>10 

20:33:30.81 

41:15:22.7 

Bllb 

MT91 

18 

556 


6.47 

5.53 

4.39 

28.9 

6.68 

20:31:56.27 

41:33:05.3 

B1.5V 

MT91 


169 


4.41 

3.31 

4.36 

7.8 

0 

20:33:05.55 

41:43:40.1 

B1.5V 

MT91 


403 


5.48 

4.14 

4.36 

10.7 

>10 

20:33:34.36 

41:18:11.6 

B1.5V 

MT91 


575 


6.73 

4.49 

4.36 

12.6 

>10 

20:33:48.88 

41:19:40.9 

B1.5V 

MT91 

71 

646 


5.11 

3.83 

4.36 

9.3 

>10 

20:32:13.75 

41:27:42.0 

B1.5V 

K07 


216 


4.28 

3.61 

4.36 

8.7 

>10 

20:34:43.51 

41:29:04.8 

B1.5III 

MT91 

30 

793 


5.73 

4.52 

4.36 

12.7 

>10 

20:32:26.50 

41:19:13.7 

B1.5III+B1V 

K07,K08 


252 


5.40 

3.50 

4.36 

8.3 

4.65 

20:32:46.74 

41:26:15.9 

B1.5III 

K07 


325 


4.98 

3.39 

4.36 

8.0 

0 

20:33:18.55 

41:15:35.4 

B2Ve 

K07 

64 

488 


8.27 

4.38 

4.34 

11.6 

>10 

20:31:22.10 

41:12:03.0 

B2V 

N08 



A30 

6.27 

4.06 

4.34 

10.0 

>10 

20:32:27.30 

40:55:18.0 

B2V 

H03 



A39 

6.05 

4.35 

4.34 

11.5 

>10 

20:30:59.43 

41:35:59.6 

B2V 

K07 


042 


4.58 

3.35 

4.34 

7.7 

0 

20:31:30.49 

41:37:15.6 

B2V 

K07 


097 


4.40 

2.94 

4.34 

6.8 

0 

20:32:03.01 

41:32:30.7 

B2Ve 

K07 


186 


4.53 

3.19 

4.34 

7.3 

0 

20:32:07.95 

41:22:00.3 

B2V 

K07 


202 


4.92 

3.23 

4.34 

7.4 

0 

20:32:14.63 

41:27:40.3 

B2V 

K07 


221 


4.67 

3.43 

4.34 

8.0 

5.15 

20:32:19.66 

41:20:39.7 

B2V 

K07 


234 


4.47 

3.48 

4.34 

8.1 

8.63 

20:32:22.15 

41:27:41.7 

B2V 

K07 


241 


4.18 

3.32 

4.34 

7.7 

0 

20:32:25.50 

41:24:51.8 

B2V 

K07 


248 


4.57 

3.47 

4.34 

8.1 

8.27 


continued on next page 


© 0000 RAS, MNRAS 000, 000-000 





22 Wright et al. 


Table Bl: Massive stars in Cygnus OB2 


RA 

(J2000.0) 

Dec 

(J2000.0) 

Spectral 

type 

Spectral 

reference 

Source number 

S58 MT91 C02 

A v 

(mag) 

log L 
( L @) 

T eff 

K 

M 

(Me) 

Age 

(Myrs) 

20:32:32.68 

41:27:04.4 

B2V 

K07 

275 


3.92 

3.17 

4.34 

7.3 

0 

20:32:37.78 

41:26:15.3 

B2V 

K07 

295 


4.25 

3.21 

4.34 

7.4 

0 

20:32:42.90 

41:20:16.4 

B2V+B3V 

K07.K12 

311 


4.95 

3.29 

4.34 

7.6 

0 

20:33:22.49 

41:22:16.9 

B2V 

K07 

513 


5.07 

3.33 

4.34 

7.7 

0 

20:33:24.78 

41:22:04.5 

B2Ve 

K07 

522 


4.90 

3.34 

4.34 

7.7 

0 

20:33:31.68 

41:21:46.1 

B2V 

K07 

561 


4.67 

3.38 

4.34 

7.8 

2.13 

20:33:47.63 

41:09:06.5 

B2V 

K07 

639 


5.94 

3.60 

4.34 

8.4 

>10 

20:33:48.83 

41:37:39.7 

B2Ve 

K07 

650 


5.70 

3.35 

4.34 

7.7 

0 

20:31:56.90 

41:31:48.0 

B2III 

K12a 

174 


4.39 

3.75 

4.34 

00 

OO 

>10 

20:32:26.10 

41:29:39.0 

B2III 

K07 

250 


3.91 

3.42 

4.34 

8.0 

4.66 

20:32:27.26 

41:21:56.2 

B2III 

K07 

255 


5.06 

3.15 

4.34 

7.2 

0 

20:32:30.72 

41:07:04.1 

B2III 

K07 

264 


3.56 

3.37 

4.34 

7.7 

1.25 

20:33:11.39 

41:17:58.9 

B2III 

K07 

441 


5.18 

3.31 

4.34 

7.6 

0 

20:33:48.40 

41:13:14.1 

B2III 

K07 

645 


6.19 

3.58 

4.34 

8.4 

>10 

20:30:50.75 

41:35:06.0 

B2II 

MT91 

021 


4.42 

3.32 

4.34 

7.7 

0 

20:32:46.45 

41:24:22.4 

B2.5V 

K07 

322 


4.85 

2.96 

4.29 

6.1 

0 

20:32:31.42 

41:30:51.4 

B2.5V 

K07 

268 


5.19 

3.30 

4.29 

6.8 

>10 

20:32:54.35 

41:15:22.1 

B3V 

K07 

37 358 


6.55 

3.52 

4.25 

7.0 

>10 

20:31:33.59 

41:36:04.3 

B3V 

K07 

106 


4.58 

2.92 

4.25 

5.6 

>10 

20:31:41.60 

41:28:20.9 

B3V 

K07 

129 


4.72 

3.06 

4.25 

5.9 

>10 

20:31:55.28 

41:35:27.8 

B3V 

K07 

164 


4.54 

2.73 

4.25 

5.3 

0 

20:31:59.82 

41:37:14.3 

B3V 

K07 

179 


5.02 

3.33 

4.25 

6.4 

>10 

20:32:06.85 

41:17:56.8 

B3V 

K07 

200 


6.41 

3.87 

4.25 

8.0 

>10 

20:32:15.03 

41:19:30.8 

B3V 

K07 

222 


5.02 

3.00 

4.25 

5.8 

>10 

20:32:38.34 

41:28:56.6 

B3V 

K07 

298 


4.47 

2.95 

4.25 

5.7 

>10 

20:33:33.38 

41:08:36.3 

B3V 

K07 

568 


6.71 

3.63 

4.25 

7.3 

>10 

20:31:34.12 

41:31:08.0 

B3IV 

K07 

108 


4.53 

2.64 

4.25 

5.1 

0 

20:32:04.74 

41:28:44.5 

B3IV 

K07 

191 


3.83 

2.94 

4.25 

5.6 

>10 

20:32:49.67 

41:25:36.4 

B3III 

K07 

336 


4.23 

2.97 

4.25 

5.7 

>10 

20:33:33.97 

41:19:38.4 

B3I 

K07 

573 


5.25 

3.58 

4.25 

7.1 

>10 

20:32:40.88 

41:14:29.3 

B3.5Ia+ 

MT91 

12 304 


10.18 

6.09 

4.14 

110 

3.0 

20:32:21.76 

41:34:24.6 

B4V 

K07 

239 


3.92 

2.63 

4.22 

4.7 

2.66 

20:32:32.34 

41:22:57.6 

B4V 

K07 

271 


5.14 

3.02 

4.22 

5.4 

>10 

20:33:10.27 

41:23:44.9 

B4III 

K07 

427 



Background 


20:31:56.23 

41:35:12.3 

B5V 

K07 

170 



Background 


20:32:32.54 

41:26:46.7 

B5V 

K07 

273 


4.53 

2.49 

4.20 

4.3 

3.92 

20:33:11.81 

41:24:05.8 

B5V 

K07 

444 


5.25 

3.22 

4.20 

5.7 

>10 

20:33:13.37 

41:26:39.7 

B5V 

K07 

453 


4.27 

2.71 

4.20 

4.7 

>10 

20:33:14.34 

41:19:33.0 

B5V 

K07 

459 


6.65 

3.45 

4.20 

6.2 

>10 

20:33:47.58 

41:29:57.7 

B5V 

K07 

641 


5.71 

3.39 

4.20 

6.1 

>10 

20:33:19.26 

41:24:44.8 

B5IV 

K07 

493 


5.56 

2.98 

4.20 

5.2 

>10 
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